arXiv:1509.06770v4 [astro-ph.CO] 6 Apr 2016 


Prepared for submission to JCAP 


Robust forecasts on fundamental 
physics from the 
foreground-obscured, 
gravitationally-lensed CMB 
polarization 


Josquin Errard,* ®’^'^’^ Stephen M. Feeney,* ® Hiranya V. Peiris,^ 
and Andrew H. Jaffe® 

“Sorbonne Universites, Institut Lagrange de Paris (ILP), 

98 bis Boulevard Arago, 75014 Paris, France 
''LPNHE, CNRS-IN2P3 and Universites Paris 6 & 7, 

4 place Jussieu, F-75252 Paris Cedex 05, France 

'^Computational Cosmology Center, Physics Division, Lawrence Berkeley National Lab, 
1 Cyclotron Road, Berkeley CA 94720, USA 
'^Space Sciences Laboratory, University of California, 

7 Gauss Way, Berkeley CA 94720, USA 

^Department of Physics, Blackett Laboratory, Imperial College London, 

Prince Consort Road, London SW7 2AZ, United Kingdom 
•^Department of Physics and Astronomy, University College London, 

Gower Street, London WCIE 6BT, United Kingdom 

*these authors contributed equally to this work 



E-mail: josquin.errard@lpnhe.in2p3.fr, s.feeney@imperial.ac.uk, 
h.peiris@ucl.ac.uk, a.jaffe@imperial.ac.uk 

Abstract. Recent results from the BICEP, Keck Array and Planck Collaborations demon¬ 
strate that Galactic foregrounds are an unavoidable obstacle in the search for evidence of 
inflationary gravitational waves in the cosmic microwave background (CMB) polarization. 
Beyond the foregrounds, the effect of lensing by intervening large-scale structure further ob¬ 
scures all but the strongest inflationary signals permitted by current data. With a plethora 
of ongoing and upcoming experiments aiming to measure these signatures, careful and self- 
consistent consideration of experiments’ foreground- and lensing-removal capabilities is 
critical in obtaining credible forecasts of their performance. We investigate the capabilities 
of instruments such as Advanced ACTPol, BICEPS and Keck Array, CLASS, EBEXlOK, 
PIPER, Simons Array, SPT-3G and SPIDER, and projects as COrE-|-, LiteBIRD-ext, 
PIXIE and Stage IV, to clean contamination due to polarized synchrotron and dust from 
raw multi-frequency data, and remove lensing from the resulting co-added CMB maps (ei¬ 
ther using iterative CMB-only techniques or through cross-correlation with external data). 
Incorporating these effects, we present forecasts for the constraining power of these exper¬ 
iments in terms of inflationary physics, the neutrino sector, and dark energy parameters. 
Made publicly available through an online interface, this tool enables the next generation of 
CMB experiments to foreground-proof their designs, optimize their frequency coverage to 
maximize scientific output, and determine where cross-experimental collaboration would be 
most benehcial. We hnd that analyzing data from ground, balloon and space instruments in 
complementary combinations can significantly improve component separation performance, 
delensing, and cosmological constraints over individual datasets. In particular, we find that 
a combination of post-2020 ground- and space-based experiments could achieve constraints 
such as cr(r) ~ 1.3 x 10“'^, cr(nt) ~ 0.03, a{ns) ~ 1.8 x 10“^, (t(q;s) ~ 1.7 x 10“^, 
a{My) ~ 31 meV, a{w) ~ 0.09, a{wo) ~ 0.25, a{wa) ~ 0.50, a{Nes) ~ 0.024 and 
(T(Dk) ~ 1.5 X 10“^, after component separation and iterative delensing. 
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1 Introduction 

Cosmic microwave background (CMB) B-mode polarization, generated by gravitational 
lensing of E-mode polarization at arcminute scales [1] and, potentially, by inflationary 
gravitational waves on degree scales [2, 3], provides a unique window into the physics 
of both the early and evolved universe. As an integrated measure of the geometry and 
growth of structure in the Universe, the measurement of lensing B modes provides both 
insight into late-time physics, such as the damping of structure formation by massive 
neutrinos [4, 5] and the influence of dark energy [6-8], and an opportunity to break the 
CMB geometric degeneracy [9, 10] and hence constrain the curvature of the Universe 
without resorting to external datasets [11]. The detection of degree-scale B modes would 
provide convincing evidence that the Universe entered an inflationary phase at very early 
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times [12], with measurement of their amplitude constraining the energy scale at which 
inflation took place [13, 14]. The lure of such scientific potential has driven considerable 
effort into the development of dedicated i?-mode observatories, and this effort has now 
begun to come to fruition. The past year has yielded the first measurements of the lensed 
i?-mode polarization [15-20], though primordial B modes remain elusive. The drive for ever 
more sensitive i?-mode measurements is therefore ongoing, with numerous ground-based, 
balloon and satellite missions funded and in various stages of taking data, and many more 
proposed.^ 

A natural consequence of the low amplitude of the i?-mode polarization is that obser¬ 
vatories designed to detect B modes can also make exquisite measurements of the scalar- 
sourced £'-mode polarization. Cosmic-variance limited E modes, and indeed their cross¬ 
correlation with temperature, provide even more constraining power than their tempera¬ 
ture counterparts [21]. The next generation of CMB polarization experiments will therefore 
tighten our constraints on the standard cosmological parameters by factors of ~ 3, as well 
as constrain additional parameters such as the number of relativistic degrees of freedom. 

However, to fully exploit the constraining power of CMB polarization one must over¬ 
come several obstacles. Recent observations [16, 22, 23] have emphatically demonstrated 
that next-generation CMB experiments will have to characterize, control and remove sig¬ 
nificant astrophysical foregrounds. Polarized Galactic foreground contamination turns out 
to be important at all Galactic latitudes, either due to large total intensity or high polar¬ 
ization fraction. Beneath the astrophysical foregrounds, the lensing B modes — though 
highly important for cosmological constraints — are a significant limiting factor in de¬ 
tecting low values of the tensor-to-scalar ratio [24-26]. Delensing the B modes [25-33] 
is therefore a further necessity for maximizing the scientihc output of the most sensitive 
polarization observatories. Significant prior effort has gone into CMB experimental design 
including the mitigation of foreground contamination and delensing [34-43]. The Planck 
Collaboration has recently released maps of polarized dust and synchrotron [23], strongly 
motivating an updated evaluation of the abilities of current and planned CMB observatories 
to address these challenges. In this paper, we employ a parametric maximum-likelihood 
approach (Ref. [44] and references therein) to estimate the ability of a range of pre- and 
post-2020 CMB experiments to remove polarized foreground contamination, delens the re¬ 
sulting B modes, and constrain the inflationary, neutrino, and dark energy sectors. The 
tool developed for this work allows the specification of instrumental configurations (ob¬ 
served frequencies, noise levels, beam sizes and sky fractions; it does not currently model 
the effect of further instrumental systematic errors) and uses a Fisher-matrix framework 
to forecast errors on the desired parameters. It has been made publicly available on a 
Lawrence Berkeley National Laboratory (LBNL) machine,^ with a user-friendly web inter¬ 
face allowing any specific instrumental configuration to be studied. 


^See, for example, http://lambda.gsfc.nasa.gov/product/expt/. 
^turkey.Ibl.gov 
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Section 2.5, Fig.5, 
Appendix A 



Figure 1. Schematic of the adopted methodology described in Sect. 2. 


The paper is organized as follows. Section 2 introduces the component separation 
approach, the delensing technique, and the Fisher framework that we use to evaluate the 
scientific performance of different experimental configurations. Section 3 summarizes our 
results, and Sect. 4 describes the web interface which is publicly released with this work. 
Finally, we discuss the implications of our results in Sect. 5. In addition. Appendix A details 
the specifications of each instrument considered in this work, and Appendix B provides a 
comprehensive set of their forecasted constraints. 


2 Methodology 

We first describe the method we use to evaluate the scientific performance of a range 
of CMB experimental setups in the presence of Galactic foregrounds and gravitational 
leasing of the CMB. The precise specifications assumed for the experiments are described 
in Appendix A. Based on these specifications, we estimate each experiment’s ability to first 
separate Galactic foreground emissions and then gravitational leasing from the primordial 
CMB, and hence constrain the cosmological parameters of interest. Figure 1 summarizes 
the methodology; in the subsections below, we describe the main steps in detail. 
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Figure 2. Left panel: Angular power spectra showing primordial B modes, leasing B modes, 
total intensity, and E modes, as well as the total contribution of polarized B-mode foregrounds 
(dust plus synchrotron), expected on the cleanest 1-90% of the sky, at 100 and 200 GHz. Note 
that, as these results are derived from Planck's Galactic masks and are not therefore optimized for 
high-resolution, ground-based instruments, there is potential for discovery of small patches of sky 
(e.g., /sky ^ 5%) cleaner than those indicated here. Right panel: The ratio of power spectra of 
foreground and lensing B modes to primordial B modes, assuming a tensor-to-scalar ratio r = 1. 
The contours indicate, in effective values of r, the contamination due to foregrounds and lensing on 
primordial B-mode measurements. The x- and y-axes correspond to the multipole i and frequency 
of observation, in GHz, respectively. The level of input foregrounds are estimated on a 50% patch 
of the sky. 


2.1 Foreground removal 

As illustrated in Fig. 2, polarized contamination from astrophysical foregrounds is an un¬ 
avoidable challenge in the quest for primordial B-mode measurements. The left panel 
shows CMB temperature, B-mode and B-mode angular spectra for different values of the 
tensor-to-scalar ratio r, as well as the expected amplitude of dust and synchrotron B 
modes, estimated using various portions of the Planck data [23]. The right panel shows, 
as a function of frequency of observation and multipole i, the ratio of foreground plus 
lensing B modes to primordial B modes, assuming a tensor-to-scalar ratio r = 1. Outside 
the half-sky Galactic mask that was considered, the minimum contamination is reached 
at multipole ~ 80 and frequency ~ 74 GHz (as confirmed by Ref. [45, 46]). This 
foreground-minimum region has an effective foreground amplitude of r ~ 0.1. Reaching 
values of r less than 0.1 (or, indeed, improving the significance of any future detection of 
r of this magnitude) therefore requires highly effective cleaning of both foregrounds and 
lensing B modes. 
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2.1.1 Formalism 

Our approach is based on the parametrization of the emission laws for two important 
astrophysical contaminants, namely dust and synchrotron [44, 47-49].^ To estimate the 
impact of performing component separation with a given instrument, we use a parametric 
maximum-likelihood component-separation approach, as implemented in Ref. [47]. We 
assume the following linear data model, in which the signal measured in each frequency 
channel i at each pixel p is given by 

di{p) = Aij Sj(p) + ni(p) , (2.1) 

where we assume summation over the repeated index j that runs over the underlying 
components (i.e., CMB, dust and synchrotron), and 

• d(p) is a multi-frequency data vector, with each entry corresponding to a different 
frequency channel; 

• s(p) is a multi-component sky signal vector, each entry of which corresponds to a 
different polarized sky component to be estimated from the data; 

• A is a mixing matrix defining how the components need to be combined to give a 
signal for each of the considered frequency channels; and 

• n(p) is a vector containing the instrumental noise in each frequency channel, which 
is assumed to be Gaussian and uncorrelated with dispersion N. 

We will switch between index (Aij) and matrix (A) notation as convenient. For simplicity, 
both A and N are assumed to be pixel independent. The mixing matrix A takes into 
account each experiment’s bandpasses 

Aij = j 5i{i')Af'^(i^,i/y.ef)du, (2.2) 

where 6i(iy) is the fractional bandpass of the frequency channel, defined such that 
f 5i(v)du = 1, and A''®'"(z^, t'ref) is a vector describing the frequency dependence of the 
components before bandpasses are taken into account. 

®We restrict ourselves here to the polarized foregrounds known to be significant at CMB frequencies. 
As measurements of the polarized CMB sky improve, they may reveal further complexity in both the form 
and number of foreground components, with potential additional sources including spinning and magnetic 
dust [50, 51] and carbon monoxide [52]. We leave the treatment of these potential foregrounds to future 
investigation. 

"'Note there exists a complementary approach to component separation that removes foregrounds by 
combining multi-frequency data to minimize the variance of the CMB component (see, e.g., Ref. [53] for the 
various methods applied to Planck data). While we do not explicitly treat such methods in this analysis, 
we can expect broadly the same performance as for parametric models. 
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We assume that component separation comprises two independent steps. First, we 
estimate the mixing matrix, A, from the various frequency maps produced by a given 
instrument. In the parametric framework assumed here, estimating the mixing matrix 
corresponds to determining the parameters describing the frequency dependence of the 
dust and synchrotron emission. Second, we linearly combine the different frequency maps 
to disentangle dust and synchrotron from the CMB — inverting Eq.(2.1) to obtain the 
sky signals s(p) — using the estimated mixing matrix. This process leads to a noise 
level in the reconstructed map which is higher than the simple quadratic combination of 
sensitivities from all frequency channels. Although the first step is specific to a considered 
component separation method, the second step is general to any approach — though it 
can be performed either in pixel or Fourier space. 

Estimation of the mixing matrix A. As supported by the latest Planck results 
[23], we assume for the synchrotron frequency dependence a simple power-law form with 
spectral index /3s, 

^ (2.3) 


where the reference frequency = 150 GHz. We consider a modified grey-body emission 
law for the dust, 

n ,1 ^^ref 

_ r u e - 1 

Vz^ref / 




(2.4) 


gfcTj _ I 

where /3d is the dust spectral index and Td the dust temperature. As the dust temperature 
and spectral index are strongly anti-correlated in the frequency ranges targeted by CMB 
experiments, with higher frequencies required to constrain the temperature, we choose to 
fix the temperature to a reasonable value and only fit for the spectral index in the following. 
The component-separation scenarios studied in this work are summarized in Table 1. We 
consider the most realistic scenario of synchrotron, dust and CMB to derive our main 
results, though in several intermediate figures we consider simplified foreground models for 
pedagogical purposes. 

Instead of performing simulations and running a full likelihood analysis, we use a 
Fisher-matrix approach to estimate the uncertainty on each foreground component’s spec¬ 
tral index. As derived in Ref. [44], the second-order derivatives of the j3 spectral log- 
likelihood, T(/3), taken at the true values of the parameters, can then be written as 



true /3 

N-^A (A'^N’^A)”^ 


5/3' 


(2.5) 


d/3 5/3' 




where s(p) are template maps for each desired component (typically CMB, dust and syn¬ 
chrotron); see Sect. 2.1.2. S is a matrix of size where Upar is the number of free 
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Table 1. Summary of the component-separation scenarios and approaches to spatially-varying 
spectral indices described in Sect. 2.1. 


label 

no comp. sep. 

d 

s 

s-bd 


CMB 

CMB 

CMB 

CMB 

components 

only 

+ dust 
only 

-|- synchrotron 
only 

+ synchrotron 
+ dust 


label 

Up approach 

A-expansion approach 

description 

Assume Upatch = 12 x /^ky x 4^ 
independent patches within 
observed sky fraction /sky. 
Assume this hypothesis 
artihcially boosts foreground 
residuals as oc ripatch- 

Extend mixing matrix A and sky 
signal s(p) by expanding emission 
laws with respect to spectral 
indices /3d and /3s around 
their mean values. Noise after 
component separation is boosted. 


parameters in the mixing matrix A. In our case, A = A(/3d,/3s) and ripar = 2. The 
uncertainty on the estimation of the foreground spectral indices is then given by [48] 

a(/3) (2.6) 

The residual foregrounds left in the CMB map after component separation can be derived 
from this uncertainty; their power spectrum is given by 

( 2 . 7 ) 

k,k' j,j' 


where are the input foreground spectra with j, j' G {cmb, dust, synchrotron}^. The 
element is as defined in [44] 


jj' _ Oj Oj' 


with = — 




JOi 


( 2 . 8 ) 

(2.9) 


The foreground residuals can ultimately bias the estimation of primordial B modes in the 
reconstructed CMB map, see Ref. [54]. To study the potential biases due to imperfect 
component separation, we therefore define the effective amplitude of foreground residuals, 

®Note that cross-spectra, and in particular (7)’“'** ^ ^j.g in the sum of Eq. (2.7). 
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1=20 


red 


^^200 ^fgres 

l^£=20 '-'£ _ 

^;p200 ^SB.prim/ _ -, x 

Z^£=20'^l 


( 2 . 10 ) 


( 2 . 11 ) 


Equation (2.5), through the term s{p) s^(p), contains the a priori assumption that 
the foreground spectral indices are constant over the sky. This is, however, not necessarily 
the case [23], and we account for this assumption in two ways. Firstly, following Ref. [49], 
we break each experiment’s observable patch into ripatch independent regions, within which 
the foreground spectral indices are assumed to be constant. Apportioning the data this 
way naively leads to the approximation (j{j3) oc y^Upatcbj or equivalently oc ripatch) 

degrading the performance of the component separation and boosting the impact of the 
foreground residuals. Following Refs. [22, 23], we take the size of these independent patches 
to match Aside = 4 HEALPix pixels [55], requiring 


^patch — [12 X /sky X 4 J 


( 2 . 12 ) 


unique spectral indices to be fitted by an experiment covering a sky fraction /sky (rounding 
down to the next lowest integer). This approach, referred to as the “up approach” in the 
following, is the method used to derive the main results of this work. 

Following Ref. [56], we also consider a second approach utilizing the first-order ex¬ 
pansion of the mixing matrix A around the mean values of each foreground parameter 


A(/3)« A(/3) + 5/3(p) 


dA 

w 


+ 0{6l3{pf), 

/3 


(2.13) 


where (3 = {(3)p is computed over the observed sky pixels p, 5j3 = j3 — j3 and, as before, the 
dust temperature Td is assumed constant across the sky. By increasing the dimensionality 
of the mixing matrix A, we can extend the parametric framework to account for spatially 
varying foreground parameters. Within this extended framework, the sky signal used in 
Eq. (2.1) is 


®(p) — [Scmb(p)) 'Sdust(p)) ^Pdip) Sdust(p) ; ^sync (p), 6f3sip) •^sync ip)] , 

where each Sj{p) G {Q,U}, and the corresponding mixing matrix A reads 

5Adust . 9As 


A = 


cmbi -^dust) 


9/3d 


/A, 


‘-sync 


sync; 


d/3s 


(2.14) 


(2.15) 


We denote this method the “A-expansion approach”. As we shall see in the following 
section, the expansion of the mixing matrix necessarily results in higher noise levels after 









component separation; however, the addition of new degrees of freedom in the foreground 
modeling can potentially yield better control of foreground residuals, provided an experi¬ 
mental configuration contains enough frequency channels to fully characterize the compo¬ 
nents present. An additional source of bias not considered here could result from using 
the wrong foreground model to clean the CMB. Throughout this work we assume that the 
foreground model is correct; see Ref. [42] for a low-multipole treatment of this systematic. 

The question of which of the Up and A-expansion approaches is most appropriate 
is difficult to answer currently, as the evidence for spatial variation of the foreground 
parameters [23] is still inconclusive. As yet, there is no compelling evidence for a particular 
angular scale in the variation of the spectral indices, and each project will have to determine 
whether the Up approach is a good approximation for their data. This will depend on the 
chosen region of the sky and on the observational integration depth. The current generation 
of (relatively) noisy small-patch experiments will likely be driven towards adopting the Up 
approach—and hence relying on priors derived from foreground observations—in order to 
retain constraining power. The increased sensitivity of future projects, coupled with their 
broader frequency and sky coverage, will make the A-expansion approach feasible, allowing 
critical tests of the priors and parametrizations used in the Up approach. 

Linear combination of frequency maps given A. Once the parameters of the 
mixing matrix are estimated, and assuming the estimates are unbiased, the noise variance 
in each pixel of the reconstructed CMB map in (/rK-arcmin)^ is given by [47] 


2 _ 2 

f^CMB = 7p 


(A^N-^ 


A) 


-1 


J CMB CMB 


(2.16) 


where 7 p is the angular size of a sky pixel in arcmin. Due to the A^N^^A inversion, cjcmb 
is necessarily larger than the quadratic combination of all the sensitivities per channel.® 
The extent of this noise degradation depends on the condition of the mixing matrix A, 
and thus the larger mixing matrix employed in the A-expansion approach, Eq. (2.13), can 
lead to worse noise levels than the Up approach. Complementary to the effective amplitude 
of residuals reg dehned in Eq. (2.11), we introduce the dimensionless noise degradation A 
dehned as 

(2A7) 

\ <7quad / 

where Uquad is the simple quadratic combination of sensitivities across all frequency chan¬ 
nels. 

Eigure 3 shows the behavior of cjcmb and A in the simple case of an experiment with 
three frequency channels — 150, 220, and 353 GHz — observing the CMB obscured only 
by dust. The sensitivities of the two low-frequency channels are allowed to vary, but the 
highest frequency has a noise level of 200 ^uK-arcmin, assumed to come from Planck. The 

®crcMB is equal to the quadratic combination of all sensitivities only in the trivial case of A being the 
identity matrix, i.e., each component only existing at a single frequency. 
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Figure 3. Left panel: the noise in the final CMB map for an experiment with two frequency 
channels of varying polarization sensitivity centered on 150 and 220 GHz, combined with Planck's 
353 GHz channel. The sole foreground contaminant is dust. Right panel: the noise degradation 
factor A between the final CMB map and the quadratic combination of all channels for the same 
experimental configuration. 


noise in the reconstructed CMB map (left panel) is limited by the 150 GHz sensitivity 
across most of the parameter space considered, and we conclude that 150 GHz sensitivity 
should therefore be prioritized over 220 GHz until uiso < 5 /rK-arcmin. Below this value, 
there is an optimal balance to be found between the sensitivities in both 150 and 220 GHz 
channels. 

Turning to the noise degradation plot (right panel), we see exactly how much sensi¬ 
tivity is wasted by investing in 220 GHz sensitivity at the expense of 150 GHz: the noise 
in the final CMB map is over 20 times larger than the raw quadratic combination in the 
top-left of the plot. The noise is degraded the least when the 220 GHz noise is roughly 
twice that at 150 GHz. This toy example shows that Planck's 353 GHz band can be a 
good dust-monitoring channel, though it is important to note that the foreground resid¬ 
uals resulting from imperfect estimation of the mixing matrix have not been considered 
here. Furthermore, adding synchrotron to the problem or letting the spectral indices in A 
vary across the sky would significantly degrade ctcmb (and therefore A) if more frequency 
channels with reasonable sensitivities were not available. 

2.1.2 Input data 

The computation of the uncertainty on the foreground spectral indices, Eq. (2.5), requires 
s(p) — polarized CMB, dust and synchrotron template maps — as well as the spectral in¬ 
dices Pdip) and /3s{p) in the case of the A-expansion approach. We use Planck component- 
separated maps [23, 53] for this purpose, degrading all maps to their highest common 
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resolution, HEALPix iVgjde = 128. Following Refs. [23, 57], we further assume 


/3d = 1.59 

rd = 19.6K (2.18) 

/3s = -3.1 


to build the elements of the input mixing matrix A. Throughout this work, we set the dust 
temperature to its true value, and vary the two spectral indices /3d and /3s. Where Planck 
is not directly included in the experimental configuration, we additionally adopt Gaussian 
priors from Planck's constraints on the foreground parameters [23, 57]. We transform 
Eq. (2.5) as 5] —)■ S) + ^pianck^ assuming a diagonal ^pianck such that 


(TPlanckiPd = 1.59) = 0.04 (2.19) 

PlanckiPs — 3.1) = 0.4. 


Regarding the foreground power spectra, as 
power law forms for dust and synchrotron, 

oc 


suggested by Ref. [22], we use the following 

(f)^'‘ 




oc — 


2.6 


where £o = 80. Several of these input spectra are shown in Fig. 2. Note that the spatial and 
angular dependencies only enter into the calculation of the residual power spectrum, 

Eq. (2.7). As we shall see, in many cases the experimental performance is dominated by the 
post-component-separation noise, ctcmb, which does not assume any information on the 
foregrounds beyond their frequency dependence.^ Finally, all the sky patches that we study 
throughout the paper use the Galactic masks created for the Planck data analysis [58].® 


2.2 Delensing 

For experiments with foreground-cleaned noise levels below 4-5 /xK-arcmin, GMB lensing 
becomes the dominant obstacle to observing weak primordial B modes [24] . Removing this 
contaminant through delensing [25-27] requires a measurement of the lensing potential, 

^Note that Planck has estimated the spectral indices in Eq. (2.20) with uncertainties of ~ 5 %, sourced by 
both genuine variation of the spectral indices across the sky and uncertainty at a given position. Varying the 
spectral indices within this range of uncertainty can—depending on the instrumental setup considered—lead 
to percent-level variations of the effective amplitude of foreground residuals, reff, and sub-percent variations 
of constraints on cosmological parameters sensitive to large angular scales, such as r. In the remainder of 
this work we will neglect this uncertainty. 

®These apodized Galactic masks and all other Planck data products can be downloaded from the Planck 
Legacy Archive at http://pla.esac.esa.int/pla/. 
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which one can use to estimate the lensed CMB B modes for subtraction from the total 
observed signal. We base our delensing procedure on Ref. [31], which provides the following 
analytical expression for the estimated leasing B modes: 


est 


1 

2£ + l 


E l fEB 

\hhi2 


X 


(Ct*f 



( 2 . 21 ) 


where 


is a geometric coupling factor. The delensed B mode is then given by 


^BB,del _ ^BB^lens ^BB,est 

~ 


( 2 . 22 ) 


Note that the presence of the noise in Eq. (2.21) guarantees that ^ 

In this work, we consider three sources for the leasing potential estimate: the CMB 
polarization itself (“CMBxCMB” delensing), the cross-correlation of the CMB and the 
cosmic infrared background (“CMB x CIB”), and measurements of the large-scale structure 
using, for example, cosmic shear or 21cm radiation (“CMBxLSS”). In the CMB x CMB 
case, the noise on this estimate is [28] 



1 


21 + 1 


Ei/f 






CBb + nBb 


i^£2 ) 

C+TNp 


-1 


(2.23) 


The crucial observation here is that the lensed B mode appears as a source of noise in the 
leasing potential estimator. Its performance therefore improves if the delensed B mode 
is substituted for the lensed B mode, in turn allowing better estimates of the leasing B 
mode, and so on. Iterating this estimator can greatly improve the ability of an experiment 
to delens the CMB [29]: indeed, in the noiseless limit perfect delensing is achievable. For 
realistic instrumental configurations, this process converges after a few steps, where we use 
the convergence criterion as below: 


E 


T\j4'4‘d 


_ ]\T4‘4'>i—£ 


< 1 %. 


(2.24) 


CMBXCIB and CMBxLSS delensing rely on non-CMB datasets for their leasing 
potential estimates, a feature that is particularly useful for experiments with moderate 
polarization sensitivities. The CIB is known to be a good tracer of the leasing poten¬ 
tial [59-62], and Refs. [32, 33] have recently demonstrated its utility in delensing the CMB. 
Following their example, to simulate CMB x CIB delensing we make the replacement 

Cf^ + Cf’/pl (2.25) 

in Eq. (2.21), where pi is the £-dependent correlation coefficient between the CMB leasing 
and CIB, as presented in Fig. I of Ref. [33]. 
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Table 2. Delensing options employed in this work and references to implementation details. 


label 

no delensing 

CMB X CMB 

CMB X CIB 

CMB X LSS 

references 

- 

Ref. [31] 

Ref. [33] 

Ref. [31] 

properties 

no delensing 
performed 

iterative delensing 
using CMB 
polarization only 

delensing using 
correlation between 
CMB and CIB 

delensing with 
perfect LSS 
measurements 

to Zmax 


The CMBxLSS delensing case requires high signal-to-noise large-scale structure mea¬ 
surements extending to high redshift (e.g., cosmic shear [30] or 21cm radiation [63-65] 
measurements). Following Ref. [31], we assume a (crude but illustrative) sharp cutoff in 
LSS tracers at Zmax, providing a perfect estimate of the gravitational potential out to this 
redshift. In this setting, the estimated leasing B mode takes the same form as Eq. (2.21), 
with the following replacement 




cf^ + 





(2.26) 


where is the power spectrum of the lensing potential including contributions from 

-2 < -Zmax Only (calculated using a modified version of CAMB [66]). In this case, delensing 
performance is limited by the noise in the E measurements and the redshift range of the LSS 
measurements. In this work, we take Zmax to be 3.5: the proposed limit of the capabilities 
of the Large Synoptic Survey Telescope (LSST) [67]. 

To evaluate the quality of the various delensing techniques in each experimental con¬ 
text we introduce the delensing factor a, adapted from Ref. [33] and analogous to Eq. (2.11): 


y—v200 

_ 2^i=20 _ 

y-^200 f^BB, lens ' 
1^1=20 '-'i 


(2.27) 


A summary of the delensing options is presented in Table 2. 

The left panel of Fig. 4 shows the delensing factor a. introduced above as a function 
of the post-component-separation sensitivity of a CMB experiment with a Gaussian beam 
with full-width at half-maximum (FWHM) of 3b Plotted are all three delensing methods 
considered in this work: iterative delensing using the CMB polarization (Eqs. (2.23) and 
(2.24)), and via cross-correlation with the GIB (Eq. (2.25)) and LSS (Eq. (2.26)). The first 
point to note is that delensing using the CMBxLSS correlation only slightly outperforms 
CMB X GIB delensing, even though CMBxCIB delensing makes use of existing Planck data 
products [33]. Moreover, for both cross-correlation methods the sensitivity of the CMB 


- 13 - 













Figure 4. Left panel: the delensing improvement factor, a (defined in Eq. (2.27)), for iterative 
CMBxCMB (dark orange), CMBxCIB (light orange) and CMBxLSS (dark red) delensing, as a 
function of the noise in the CMB map after component separation. We assume FWHM=3' and 
a multipole range of ^2n/fs]^y < £ < 3000. Figure 3 and its accompanying discussion describe 
three-frequency experimental configurations producing comparable noise levels after the removal of 
dust. Right panel: One-sigma limit on r = 0, cr(r = 0), as a function of both the noise in the CMB 
map and the observed fraction of the sky. Grey (black) contours correspond to the no delensing 
(CMB X CMB delensing) case. No foreground residuals are taken into account in this result. 


instrument does not matter in the regime <5 — 10 /iK-arcmin, as the £^-mode observations 
are largely cosmic-variance limited for these sensitivities (see Eq. (2.21)). Iterative delens¬ 
ing using CMB polarization becomes more effective than the cross-correlation techniques 
when the post-component separation noise level drops below ~ 3 /rK-arcmin.® This is a 
consequence of Eq. (2.23), which continues to decrease as the instrumental sensitivity 
increases, in the framework of iterative delensing [29]. 

One final point is to note that the curves become thicker as the noise level decreases. 
This effect is due to varying the area of sky to which the experiment has access (sky 
fractions considered here span the range 1-100%). Larger-scale experiments have access to 
more multipoles (we assume £rmn oc ) which allows for improved iterative delensing 
at low noise levels. 

2.3 Scientific performance forecast 

We adopt a Eisher matrix approach to provide estimates of the scientific potential of each 
experimental configuration [68]. Following, e.g., Refs. [69, 70], we write the Fisher matrix 

® Other quadratic combinations of CMB fields — and, indeed, their minimum-variance combination [27] 
— will outperform the “EBEB’’’ estimator considered here in the high-noise regime. 
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element Fij for CMB spectra as 


F- ■ — 


E 


2i+l 


/skytr 


^ de, 


(2.28) 


where 9i and 6j are two of the cosmological parameters of interest. The covariance matrix 
C, is defined as 


Cp + xp 

cr 

0 

cr 1 

f^TE 

Cf^ + 

0 

f-iEd 

0 

0 

cf ^ + X^B 

0 

nTd 

f-tEd 

0 

cf + xf_ 


(2.29) 


where Ci are the various auto- and cross-power spectra of the CMB temperature (T), 
polarization (E,B) and deflection (d) components. So as to not double-count the lens- 
ing information encapsulated in the deflection field, we use only unlensed T, E and B 
information, as denoted by barred Cis [5]. In Eq. (2.29), the diagonal elements of the 
covariance matrix contain “catch-all” Gaussian noise terms ■ For the components 
X = {T,E,B}, this noise power spectrum accounts for the effects of instrumental noise, 
imperfect foreground removal and, in the case X = B, delensing: 


_ jyBB,inst ^fgres ^BS,del. 


(2.30) 


is defined in Eq. (2.7) and in Eq. (2.22). Interpreting the foreground resid¬ 

uals as a bias with unknown shape, we parametrize 


= A 


fg res 


X C: 


'fg res 


^0 


X 



(2.31) 


and derive all Fisher constraints on cosmological parameters after marginalizing over Afg res 
and 6fg res (when foregrounds are present, cf. Table 3). The choice of parametrization in 



— 1/2 

where ^ is the instrumental white noise level of a given frequency channel u in /tKcmb- 
rad and dFWHM,^ is the full-width at half-maximum beam size in radians. We assume fully 
polarized detectors, such that . Eq. (2.32) is only valid in its 

given format in the case of no component separation, i.e., the CMB-only case of Table 1. 
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For the realistic cases in which component separation is performed, we rescale 

using A, the noise degradation due to component separation defined in Eq. (2.17), such 

that 


N. 


XX, inst 


AA, 


XX, inst 


(2.33) 


The uncertainty on the reconstruction of the deflection field, N^'^, depends on the delensing 
procedure adopted. When the lensing potential estimate is derived solely using the CMB 
(the “no delensing” and “CMBxCMB” delensing cases), the noise is given by Eq. (2.23), 
appropriately reweighted in i. When external data are used, can be derived from the 
form of the correlations between lensing tracers and the potential (Eqs. (2.25) and (2.26)). 

The Eisher matrix (Eq. (2.28)) requires partial derivatives of the covariance matrix 
with respect to the parameters of interest. To calculate these derivatives numerically, we 
evaluate four additional sets of power spectra for each parameter. Eor parameters with non¬ 
zero values in the fiducial cosmology, the additional power spectra are calculated at 90%, 
95%, 105% and 110% of the parameter’s fiducial value; for parameters whose fiducial values 
are zero, we calculate the additional power spectra at ±0.1 and ±0.05.^® We interpolate 
between these and the fiducial power spectra at each i, using a cubic spline to get a 
smooth curve, before finally forming a two-sided numerical derivative with very small step 
size (10“^) from the interpolated power spectra. The step sizes used during the first stage 
are tuned so that derivative computations converge for all multipoles and numerical noise 
is minimal. 

The Eisher formalism allows forecasting of uncertainties either conditional on the 
other parameters taking their fiducial values or marginalized over the parameters taking 
any value. Conditional errors are given simply by the inverse of individual entries in 
the Fisher matrix, 1 / ^/Fij ; marginal errors, which we employ throughout, are given by 
inverting the Fisher matrix: 

cji = a{pi) = \/[F-i]ii. (2.34) 


As an illustration of the method described above, the right panel of Fig. 4 plots the 
one-sigma limits on r = 0, a{r = 0), as a function of the 3^-FWHM example experiment’s 
post-component-separation noise level and fraction of sky covered. Grey filled contours 
show the limits without delensing, black contours the limits with iterative CMB x CMB 
delensing. We assume only i?-mode information is used to derive the limits on r, and 
all other cosmological parameters are fixed; in this setting, the Fisher formalism reduces 


to [73] 


cj(r = 0) = 


E 


{2i ± l)/sky ( 


BS,prim/ _ 


(r = 1) 




- 1/2 


(2.35) 


exceptions to these rules are the energy density in curvature, and the sum of the neutrino 
masses, M,^. For the curvature we calculate power spectra at flk = ±0.002 and ±0.001; for Mi, we use 80%, 
90%, 110% and 120% of the fiducial value. 
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where the total i?-mode noise, , comprises the post-component-separation instrumen¬ 
tal noise and the lensed B mode or its post-delensing residual, as appropriate. With¬ 
out delensing, the limits on r achievable by low-noise experiments (with noise levels < 3 
;uK-arcmin) are fixed by the sky fraction, which sets the number of modes available for 
analysis.Iterative CMB delensing completely removes this /sky performance floor, al¬ 
lowing arbitrary improvement in the limits on r = 0 as the experimental noise approaches 
zero [29]. The impact of delensing is greatest for experiments with smaller sky fractions 
as the variance due to lensing B modes becomes increasingly dominant over the putative 
signal as i increases (cf. Eq. (2.35) and Fig. 2). 

2.4 Models and fiducial cosmologies 

The aim of this work is to estimate the ability of upcoming experiments to constrain the 
parameters of extensions to the standard cosmological model, ACDM, in the inflationary, 
neutrino and dark energy sectors. Our analysis therefore always includes the standard six 
ACDM parameters: the spectral index, Ug, and amplitude, Ag, of the power spectrum of 
scalar perturbations, the optical depth to reionization, r, the Hubble constant, Hq = lOOh 
km/s/Mpc, and the current energy densities in baryons, and cold dark matter, 

The extended models are described by the running of the scalar spectral index ag, the ratio 
of the primordial tensor and scalar power spectra r, the spectral index of the tensor power 
spectrum, nt, the current energy density in curvature, Hk, the sum of the neutrino masses, 
Mjy, the number of relativistic degrees of freedom in the early Universe, A^eff; and two 
parameterizations describing the dark energy equation of state: a constant equation of 
state parameter, w, and a time-dependent version, w = w{a) = rco + (1 — a)wa, where a 
is the scale factor [74, 75]. When deriving constraints on the additional parameters, we do 
so in the context of minimally extended models: for example, constraints on r are derived 
assuming the true cosmological model can be described by the six ACDM parameters 
and r only; constraints on Aeff are derived assuming a six-parameter ACDM cosmology 
with a variable number of relativistic degrees of freedom, etc.^^ Precise definitions of the 
models considered can be found in Table 3, along with the parameter constraints for which 
the models are used. In practice, we calculate the Fisher matrix once for all parameters 
and derive constraints on individual parameters of interest by inverting the submatrices 
corresponding to the appropriate models. 

Our fiducial cosmology is a flat ACDM universe with one massive and two massless 
neutrinos. We set the majority of our parameter values to the best-fit ACDM values 
found in Planck’s analysis of a compilation of CMB, lensing, baryon acoustic oscillation, 
supernova and expansion data [76]. Specifically, we use Ug = 0.9667, Ag = 2.142 x 10“®, 

^^The stepped appearance of the contours for large sky fractions is due to the conversion of the minimum 
multipole to an integer for the sum over modes. 

^^The one exception to this rule are the constraints on rig and Oa, which are derived from a model with 
“full” inflationary freedom: a running scalar spectral index, tensors with unknown amplitude and spectral 
index, and curvature. 
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Table 3. Models considered and their parameterizations. Constraints reported on additional 
cosmological parameters are derived using specific cosmologies, as given in the right column. In 
addition, all constraints are obtained after marginalization over A{g res and big res, the foreground 
residuals’ amplitude and ^-dependence. 


r = 0.066, Ho = 67.74 km/s/Mpc, = 0.0223, = 0.1188, = Ilk = 0.0, 

Mr, = 0.06 eV, iVeff = 3.046, ij; = t(;o = — 1 and Wa = 0. For the majority of models we 
additionally assume a very small tensor component (satisfying the inflationary consistency 
relation [77], rit = —r/8), setting r = 10“^, ret = —1.25x10“'^. In order to derive interesting 
constraints on ret (using the ACDM-|-r-|-ret model) we boost this tensor component to 
r = 10“^, ret = —1.25 x 10“^. In all cases, the scalar and tensor spectral indices are 
measured at pivot scales of 0.05 and 0.002 Mpc“^, respectively. 

2.5 Instruments 

The full scope of possible instrumental configurations that could be tested by our framework 
is too large to be considered in this work, and the code is accordingly made available via a 
web interface, as described in Sect. 4. To illustrate the different options and characteristics 
of our algorithm, we focus on various experiments already in place or expected to deploy 
in two timeframes, 2015-2020 and 2020-2030, namely: 

• 2015—2020. The ground-based Stage-Ill experiments Advanced ACTPoh^, BICEP3 
and the Keck Array^^, CLASS^^, SPT-3G^®, and the Simons Array^^, as well as the 


^^Ref. [78] and E. Calabrese (private comm.). Note that an update to the Advanced ACTPol specifica¬ 
tions [79] appeared after the completion of this work. The updated design includes additional channels at 
28 and 41 GHz, which will certainly help Advanced ACTPol disentangle synchrotron from the CMB. 
Ahmed (private comm.). 

^®Ref. [80] and T. Essinger-Hileman (private comm.). 

^®Ref. [81] and B. Benson (private comm.). 

^^Ref. [82] and A. Lee and K. Arnold (private comm.). 
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Figure 5. Illustration of the frequency and multipole coverage of the instruments detailed in 
Sect. 2.5 and Appendix A. The £niax values plotted indicate the multipole at which the noise and 
beam dominate all cosmological signal (where this is lower than the maximum multipole considered 
in this work, 4000). 


balloon-borne missions EBEXIOK^®, PIPER^® and SPIDER^*^. We consider these 
experiments both alone and in combination with Planck'^^. The balloon observatories 
are designed at least partially as dust-monitoring observatories through the use of 
frequencies well above 150 GHz; Ref. [16] has unequivocally demonstrated that the 
presence of dust-monitoring channels with similar sensitivity to the CMB channels is 
crucial to the search for primordial B modes at any Galactic latitude. 

• 2020—2030. A ground-based Stage-IV-like instruments^ and the satellites LiteBIRD- 
ext^^, GOrE-|-S'^ and PIXIE^^. The ground-based Stage-IV-like instrument is assumed 
to have increased sky and frequency coverage compared to its Stage-III counterparts; 
the satellites provide high-sensitivity, broad-frequency coverage of the entire sky. We 

^®Proposal to NASA in 2015. and S. Hanany (private comm.). 

^®Ref. [83] and J. Lazear (private comm.). 

^°Ref. [84] and C. Contaldi (private comm.). 

^^Ref. [23]. 

^^In the absence of a finalized design for Stage-IV, we choose its specifications so that its noise after 
component separation, acMB, is around 1 pK-arcmin, Ref [70, 85]. 

^^http: //ltdl6. grenoble . cnrs .f r/IMG/UserFiles/Images/09_TMatsuiiiura_20150720_LTD_vl8 .pdf 
^''http: //conservancy .umn. edu/handle/11299/169642 
^®Ref. [86] and A. Kogut (private comm.). 
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forecast the performance of these instruments both in isolation and in combination. 

Figure 5 illustrates the frequency coverage and multipole range observed by each in¬ 
strument (reproducing the format used in Ref. [40]). We define the minimum observable 
multipole for each experiment using its sky patch to be t'min = [vr/(2/g,(y)], rounding up 
to the next integer. We further assume that systematic effects such as atmospheric con¬ 
tamination [87] and E-B leakage [88, 89] will render the lowest multipoles of sub-orbital 
experiments unsuitable for science, and set an floor of 20 for all such experiments (see 
Ref. [40, 42] for more optimistic forecasts). We assume a maximum multipole £max = 4000 
for all experiments, though the noise and beam for many of those considered becomes dom¬ 
inant at much lower multipoles. In Fig. 5, we plot an effective maximum multipole, chosen 
to be the lower of 4000 and the multipole at which the instrumental noise reaches 10^ 
where it dominates over even the temperature power spectrum. Note that the formalism 
does not presently account for any further systematic errors beyond those encoded in the 
choice of these parameters. 

In addition to these experiments, we also consider C-BASS^® and QUIJOTE-CMB^^ 
as dedicated synchrotron monitors when examining the more challenging A-expansion ap¬ 
proach to handling spatially varying foreground spectral indices. Due to QUIJOTE-CMB’s 
location in the Northern Hemisphere, we include it only in conjunction with space-based 
missions guaranteeing considerable overlap in coverage. The specifications for all exper¬ 
iments can be found in Tables 4, 5 and 6 of Appendix A. Note that a number of these 
experiments are still at the proposal stage, and their precise specifications are therefore 
subject to change. 

When considering the combination of two instruments A and B, to avoid double¬ 
counting data we add three Eisher matrices, the first containing information from both 
instruments from the patch of sky they have in common, the second containing any re¬ 
maining information from the instrument covering the larger sky fraction, and the third 
containing any low-I" information from the common sky patch that is only accessible to one 
experiment. The combined Eisher matrix is then 

= nfif) + F'(A/,ky) + F"io„_,(4ty (2-36) 

where /^y ^ is the common sky fraction and 

A/sky = max (/^ty, /sky) - /sty® (2-37) 

is the remainder (should there be any). To build F(/gty ®)) take the overlap in multipole 
and sky coverage of the two instruments, i.e., 

= max (2.38) 

/sky® = mia (/sky> /sky) ; 

^®http://www.astro.caltech.edu/cbass 
2'^Ref. [90]. 
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as stated previously, we consider -fmax = 4000 in all cases. With these equations satisfied, 
the frequency channels, sensitivities and beams of the two instruments are concatenated and 
fed into the foreground-cleaning, delensing and Fisher algorithms. F'(A/gky), if needed, 
is computed using the specifications of the instrument having the largest sky coverage, 
substituting A/sky for the instrument’s original sky fraction. Similarly, F"iow-£(/s^y^)) 
if needed, is computed using the specifications of the instrument observing the lowest 
multipoles, using only the multipoles between min -^min) 

As stated in Sect. 2.2, we consider two further sources of experimental data for de¬ 
lensing purposes: Planck’s 545 GHz data [33] are used in the CMBxCIB delensing forecast 
(Eq. (2.25)), and LSST-like LSS observations out to redshift Zmax = 3.5 [67] are used in 
the CMBxLSS delensing forecast (Eq. (2.26)). 

3 Results 

3.1 Foreground cleaning 

In this section, we will first illustrate the detailed foreground cleaning performance of se¬ 
lected experiments from each timeframe considered, before providing comprehensive clean¬ 
ing results for the full set of experiments. Figure 6 demonstrates the ability of two pre-2020 
and two post-2020 experiments (and their cross-correlations) to clean a variety of fore¬ 
ground components, under the np-independent patches approach. Overlaid upon the power 
spectra of primordial B modes with r = {10“^, 10“^, 10“^} (solid dark grey) and lensing B 
modes (solid light grey), we plot the post-component-separation residuals (short-dashed) 
and CMB map noise (long-dashed) for the dust-only (dark orange), synchrotron-only (light 
orange) and synchrotron plus dust (dark red) cases. The raw quadratic combination of the 
sensitivities in all channels is plotted as a dotted black curve. 

The figure is split into two panels, the top panel focusing on the pre-2020 Simons Array 
and EBEXlOK instruments, the bottom on the post-2020 Stage-IV and LiteBIRD-ext 
observatories. For the pre-2020 experiments we find noise degradation values in the ranges 

3.2 - 19, 1.6-44 and 1.6 - 5.2 for the Simons Array, EBEXlOK and their combination, 

respectively, with the largest values corresponding to the synchrotron plus dust case. For 
these experiments, the foreground residuals are subdominant to the noise, with values in 
the range 3.8 x 10“® < < 10“^, comparable in the worst-case scenario to the smallest 

cosmological signal plotted. EBEXlOK’s design as a dust observatory, amply demonstrated 
by the similarity between the pre- and post-dust-removal noise curves (A = 1.6), means it 
alone does not perform well in the cases with both dust and synchrotron. While Simons 
Array is not exceptional under any foreground assumption, its combination with EBEXlOK 
is much more robust and sensitive, though the combination is still noise- rather than 
lensing-limited for measuring primordial B modes. The increased frequency coverage and 
sensitivity of the post-2020 missions brings order of magnitude improvements in post- 
component-separation noise (1.2 < A < 2.5), with LiteBIRD-ext and PIXIE particularly 
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Simons Array x 




Figure 6. Component-separation performance for selected experiments and their cross-correlation. 
Short-dashed colored curves correspond to the residuals (Eq. (2.7)) and long-dashed colored curves 
correspond to the noise (cf. Eq. (2.16)) after removing only synchrotron (light orange), only dust 
(dark orange), or synchrotron-|-dust (dark red). These curves are compared to the primordial B- 
mode spectra for r = {10“^, 10“^, 10“^} (solid dark grey), and lensing B modes (solid light grey). 
The long-dashed black noise curve corresponds to the simple quadratic combination of sensitivities 
in all channels. Note that Planck data are not explicitly included in any of the instrumental 
configurations above; however, priors on /3(j and /3s, Eqs (2.19), derived from Planck data, are 
employed. 


robust to the foreground composition. Though the foreground residual level is at worst 
1.1 X 10“^^, the residuals become non-negligible at the largest scales attainable by these 
experiments, and these scales should therefore be treated with care. In all cases, the post- 
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component-separation noise is below the lensing, and we should therefore expect delensing 
to help considerably in constraining primordial B modes. 

Having considered a handful of experiments in detail, we now present results for the 
full suite of experiments. Figure 7 illustrates the foreground-cleaning performance of all 
experiments considered in terms of the noise degradation (A) and effective level of fore¬ 
ground residuals (reg), adopting the Up approach to spatially-varying foreground spectral 
indices. Recall that A, which is obtained by assuming the true mixing matrix A, quantifies 
the increase in variance imparted by the foreground cleaning, and the complementary rgg 
measures the level of bias left in the H-mode power spectrum due to the mixing matrix’s 
misestimation. For a given instrumental configuration, both quantities do not necessarily 
perform equally well, especially when there are too few frequency channels available. 

Concentrating first on the pre-2020 experiments. Fig. 7 indicates that, when consid¬ 
ering a single foreground component, foreground residuals are under control in most cases, 
with effective amplitudes rgg < 10“^. Similarly, the noise degradation is below 5 for most 
configurations considered. However, both the noise degradation and the residuals become 
considerably worse when including more sky components, and only the cross-correlations 
between Stage-HI instruments and Planck remain below A < 5 in the presence of dust 
and synchrotron (top-right panel). The precise values of A and rgg for the synchrotron 
plus dust case are tabulated for the cross-correlations with Planck in Appendix B. Note 
that, even when considering both dust and synchrotron, dust-monitoring instruments such 
as EBEXlOK and PIPER perform especially well in cross-correlation with Planck. As 
expected, the post-2020 instruments (bottom panel) perform better with respect to both 
criteria than the previous generation’s observatories, with their cross-correlations yield¬ 
ing even more control over foregrounds. The experiments perform consistently well, with 
typical residual contaminations of 0(10“^) for individual experiments and for 

combinations; we find noise degradation factors ranging from 2.5 for Stage-IV alone to 1.2 
for PIXIE alone and Stage-IV’s cross-correlations with PIXIE and LiteBIRD-ext. 

Eigure 8 shows the component-separation performance obtained using the expansion 
of the A matrix, as opposed to the Up approach. In this case, the foreground residuals 
remain under control (for experiments with a number of frequency channels larger than the 
dimension of the expanded A), but the noise degradation worsens considerably: we find 
noise degradation values of A > 3 even for the cross-correlations of post-2020 experiments. 
These effects arise from the added complexity of the mixing matrix: adding extra dimen¬ 
sions to A means there are sufficient degrees of freedom to fit the foreground emission 
(yielding good control of residuals), but this also leads to a larger noise boost after com¬ 
ponent separation (due to the mixing inversion in the computation of cjcmb, Eq. (2.16)). 
In order to help control the additional foreground complexity, we add C-BASS to each ex¬ 
perimental configuration as a dedicated synchrotron monitor. These results are plotted in 
Fig. 9. Adding C-BASS yields factors of at least 2-3 improvement in the noise degradation 
and almost an order of magnitude in foreground residuals for the post-2020 experiments. 
Dedicated synchrotron monitoring is necessary to maximize the scientific output of Stage- 
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Figure 7. Top panel: Component-separation performance of selected pre-2020 experiments in 
the presence of CMB, dust and synchrotron. The x-axis corresponds to the effective amplitude 
of residual astrophysical foregrounds, Teg (Eq. (2.11)), in the reconstructed CMB map. The y- 
axis is the degradation factor (Eq. (2.17)), or noise boost, driven by the linear combination of 
frequency maps during the foreground cleaning process. Upivard triangles are for satellites, circles 
are for balloons, do'wn'ward triangles are for ground-based instruments, and crosses indicate cross¬ 
correlations between experiments. Bottom panel: Same as top panel but for experiments due to 
begin their observations post-2020. Note different axis ranges. 


IV-era CMB experiments if we are to fit for spatially-varying foreground spectral indices 
(a finding confirmed by Ref. [45], who fit for scale-dependent spectral indices). We note 
that the high-frequency coverage of PIXIE is particularly effective in combination with the 
synchrotron monitoring, with PIXIE and its cross-correlations suffering noise degradations 
around half those of the other post-2020 experiments. 

As briefly mentioned in section 2.1, selecting the optimal complexity for the foreground 
parametrization amounts to trading off the bias of foreground residuals against the variance 
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A-expansion approach, post-2020 
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Figure 8. Same as Fig. 7, but using the expansion of the mixing matrix A, as described in 
Eq. (2.13). On one hand, this expansion leads to a significant increase of the noise degradation A 
due to the inversion of a larger matrix, cf. Eq. (2.16). On the other hand, the foreground residuals 
can be slightly better than for the Up approach. This is due to the modeling of both dust and 
synchrotron emission encapsulating more degrees of freedom, trading increased noise variance in 
the final CMB map for reduced residual bias. 


in the reconstructed CMB maps. Though using the A-expansion approach assumes less 
about the spatial variations of spectral indices, helping to minimize the bias, this inevitably 
leads to increased CMB noise after cleaning. Although combining data from dedicated 
foreground monitors will help reduce both bias and variance, determining whether the 
foregrounds are complex enough to warrant paying the price of constraining power will be 
difficult to check in practice, requiring careful testing using reliable sky models. 
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A-expansion approach, pre-2020 + C-BASS 
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A-expansion approach, post-2020 -I- C-BASS 
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Figure 9. Same as Fig. 8, but including C-BASS in each experimental combination. 


3.2 Delensing and detectability limits on tensor-to-scalar ratio 

As with the foreground-cleaning performance, in this section we will first illustrate the 
detailed delensing performance using our two pedagogical pre- and post-2020 experiments, 
before plotting the comprehensive delensing performance for the full set of experiments. 
Figure 10 demonstrates the ability of each experimental configuration to delens their 
foreground-cleaned CMB maps. Plotted is the fraction of the lensing B mode remaining 
after delensing, as a function of f, for CMBxCMB (solid), CMBxCIB 

(long-dashed) and CMBxLSS (short-dashed) delensing after cleaning only dust (dark or¬ 
ange), only synchrotron (light orange) or dust and synchrotron (dark red), as well as 
without performing any cleaning (black). 

From these plots, it is clear that CMB x CMB delensing for experiments in the 2015- 
2020 period (and, indeed, LiteBIRD-ext) can reduce the lensing contamination by factors 
of only a few tens of percent, with considerable variation in performance due to the compo- 
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Figure 10. Delensing performance of the different experiments. The curves show the fractional de- 
lensing residual after CMBxCMB (solid), CMBxCIB (long-dashed) and CMBxLSS (short-dashed) 
delensing. We consider cases in which, prior to delensing, synchrotron (light orange), dust (dark 
orange), or synchrotron and dust (dark red) have been removed, as well as the CMB-only case in 
which no foregrounds are present (black). Please note the different y-axis ranges. Note that Planck 
data are not explicitly included in any of the instrumental configurations above; however, priors on 
(3d and (3s, Eqs (2.19), derived from Planck data, are employed. 


sition of the foregrounds (note, however, that we only consider polarization delensing here: 
further improvement may well be possible using estimators based on the temperature or the 
minimum-variance combination of temperature and polarization [27]). Delensing through 
cross-correlation with the CIB or LSS — which require only well-measured CMB E modes 
— look promising, yielding a factor of two reduction in lensing contamination, robust 
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to the precise form of the foregrounds. Recall, however, that the B modes will still be 
noise-dominated after delensing (cf. Fig. 6), and thus we should expect only small im¬ 
provement in the attainable constraints on r = 0. Iterative CMBxCMB delensing is much 
more promising for the Stage-IV instrument, offering a factor of hve reduction in lensing 
contamination, again, robust to the foreground composition. 

Figure 11 summarizes the delensing performance of the full suite of experimental 
configurations, plotting the delensing factor a against the achievable one-sigma limits on 
r = 0 after CMBxCMB (left), CMBxCIB (center) and CMBxLSS (right) delensing. For 
clarity, only the data the two experiments have in common (i.e., the overlapping /sky and 
multipoles) is used to generate these plots. The plot is split into upper and lower panels, 
focussing on the pre- and post-2020 experiments respectively, and in all cases both dust and 
synchrotron have been removed (adopting the Up approach to spatially varying foreground 
spectral indices) prior to delensing. There are several observations to make. First, the 
limited impact of delensing on the pre-2020 experiments is clear: though around 40% more 
of the lensing can be removed by, e.g., CMBxCIB delensing than by iterative CMBxCMB 
delensing, there is minimal improvement on the one-sigma limits on r. This is, as has 
been previously mentioned, because the post-component-separation B modes are noise- 
rather than lensing-dominated. Regardless, and in particular in combination with Planck, 
the pre-2020 experiments are able to push the limits on r = 0 down to ~ 3 X 10-3. The 
post-2020 experiments’ increased sensitivity and frequency coverage, and resulting ability 
to delens, yield an order of magnitude improvement on the constraints one can place on 
r = 0; indeed, the combinations of post-2020 experiments involving a space mission will be 
able to constrain r = 0 with one-sigma errors of 10“^. 

We also study the delensing performance and detectability limits on r = 0 when 
adopting the more ambitious A-expansion approach to fitting spatially varying foreground 
spectral indices. Examining the pre-2020 experiments (not plotted), we hnd that the 
presence of Planck and C-BASS is essential in any experimental conhguration to improve 
the current limits of a{r = 0) ~ 0.1, and that the addition of further experiments improves 
the Planck plus C-BASS limits of ~ 4 x 10“^ by only a few tens of percent. The post-2020 
results, plotted in Fig. 12, are more promising. When we assume the n^-patch approach, 
the a{r = 0) limits for all post-2020 configurations are clustered in the range 1-2x10-^. 
When we constrain the spatially varying spectral indices, this range increases markedly, 
with Stage-IV in particular struggling {a{r = 0) 10 ^). The inclusion of C-BASS 

helps control the chaos, however: Stage-IV constraints improve by more than an order of 
magnitude, and all other combinations by a factor of at least two. In summary, with the 
addition of C-BASS, limits on r = 0 of 0(10“^) are achievable with post-2020 experiments, 
even with the most complex foreground composition considered. 

3.3 Full cosmological parameter constraints 

For the interested reader, a comprehensive collection of constraints on the parameters ex¬ 
tending ACDM is presented in Appendix B. Constraints are tabulated for all pre-2020 
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Figure 11. Top panel: One-sigma limits on r for the pre-2020 experiments after iterative 
CMBxCMB (left), CMBxCIB (center) and CMBxLSS (right) delensing. In all cases, dust and 
synchrotron have been removed — assuming hxed spectral indices within Up independent patches 
of sky — prior to delensing. Bottom panel: as above, but for the post-2020 experiments. 


experiments in combination with Planck, as well as for all post-2020 experiments and their 
cross-correlations. When considering Stage IV in isolation, we add Planck to provide low-£ 
information. Within each experimental table the constraints are grouped by the fore¬ 
grounds removed (none or synchrotron and dust - the Up approach is adopted by default) 
and the delensing method employed (none, CMBxCMB, CMBxCIB or CMBxLSS). For 
the post-2020 experiments, constraints are also provided for the A-expansion approach to 
foreground removal, additionally combining the experiment(s) with C-BASS and, for the 
full-sky satellites, the Northern hemisphere-based QUIJOTE-CMB. As previously stated, 
the constraints are calculated using the minimal model containing the particular parameter 
of interest, as detailed in Table 3. 
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A-expansion approach, post-2020 
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A-expansion approach, post-2020 -I- C-BASS 
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Figure 12. Top panel: One-sigma limits on r for the post-2020 experiments after iterative 
CMBxCMB (left), CMBxCIB (center) and CMBxLSS (right) delensing. In all cases, dust and 
synchrotron have been removed — accounting for spatially varying spectral indices to first order 
— prior to delensing. Bottom panel: as above, but with the addition of C-BASS as a dedicated 
synchrotron monitor. 


In this section we highlight some particularly interesting findings selected from this 
full set of forecasts. Concentrating first on the pre-2020 instruments (in combination with 
Planck), the noise degradation due to component separation is, in the majority of cases, 
less than 6, with the lowest value (A = 1.2) reached by PIPER x Planck. Foreground 
residuals are typically at the same level as the fiducial primordial B modes, reg ~ 10“^, 
though Teff < 10“® is reached by EBEXlOK x Planck and BICEPS x Keck Array x 
Planck. Though all experiments are able to remove more than 50% of the lensing signal 
by cross-correlating the CMB with tracers of the lensing potential, the constraints on r 
do not improve significantly: they are limited by the noise and residuals after component 
separation, not by the lensing. None of the pre-2020 missions can measure primordial 
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tensors with r = 0.001 — indeed, doing so would be beyond their reach even without 
foregrounds to clean — but AdvACTPol, CLASS, EBEXlOK, Simons Array and SPT-3G 
(in combination with Planck) are able to rule out r > 10“^ at the two sigma level. All 
experiments are able to place one-sigma errors on nt of ~ 0.2 (assuming r = 0.1). Though 
these constraints are not tight enough to test the inflationary consistency relation [91, 92] — 
assuming, of course, that tensors are detected at the appropriate level — they will be able 
to test alternatives to inflation such as the ekpyrotic model [93] and some models of loop 
quantum cosmology [94] at the 10-sigma level, as well as beginning to distinguish between 
standard slow-roll inflation and more general single-field models (see, e.g.. Ref [95]). When 
allowing the scalar spectral index to run, we find (T(ns) ~ 2 — 3 x 10“^, leading to a 
10-sigma measurement of the deviation from scale invariance, and cr(as) ~ 3 — 5 x 10“^, 
allowing an important test of slow roll inflation [96]. Constraints on spatial curvature of 
o'(^k) < 2.5 X 10“^ are achievable by AdvACTPol and Simons Array (in combination with 
Planck)?^ 

While delensing does not help in the search for primordial B modes, the improved 
measurement of the lensing deflection power spectrum enabled by the cross-correlation can 
significantly improve constraints on the sum of neutrino masses: AdvACTPol x Planck 
and Simons Array x Planck, for example, are able to place one-sigma limits of a{My) = 
55 meV by exploiting Planck's CIB measurements, near the threshold for distinguishing 
the degenerate from the normal and inverted neutrino mass hierarchies [85]. We stress 
that these constraints are derived from CMB (and CIB) information alone: combination 
with data from galaxy surveys such as eBOSS [97] or DESI [98] will yield even greated 
discriminatory power [99, 100]. Considering the number of neutrino species (or indeed other 
relativistic degrees of freedom), a{Nf,s) approaches ~ 0.06, a factor of ~ 6 improvement on 
current CMB-only limits [76] . In the dark energy sector, assuming a redshift-independent 
equation of state yields errors on the equation of state parameter a{w) of at least ~ 
0.16. Expanding the model to include redshift dependence understandably weakens the 
constraints, with optimal errors of ~ 0.4 and ~ 0.7 on wq and Wa, respectively. We would 
expect similar error-bar expansion for other models of dark energy with a time-varying 
equation of state. 

Turning to the post-2020 instruments, the noise degradation due to component sepa¬ 
ration, assuming the Up approach, is less than 2.5 in all cases, and the foreground residuals 
at least an order of magnitude beneath the fiducial tensor signal (though not necessarily 
the noise: see Fig. 6); these experiments efficiently and thoroughly clean polarized fore¬ 
grounds. Post-component-separation noise levels are low enough for delensing to make a 
clear impact in all cases, with any combination involving Stage-IV or COrE-|- in particular 
benefitting from CMB x CMB delensing. In the inflationary sector, one-sigma error bars on 
r = 10“^ of 1.8 X 10“^ (a 5-sigma measurement) are achievable with Stage-IV x Planck, 

^®Note that we restrict ourselves to considering CMB x CMB and CMB x CIB delensing when reporting 
the pre-2020 results here, as the observations required for CMBxLSS delensing will not be available on this 
timescale. Results for all delensing options are tabulated in Appendix B. 


- 31 - 



dropping to 1.3 — 1.4 x 10“^ in combination with any satellite mission. Testing the single- 
field inflationary consistency relation is out of reach even with r = 0.1, with the minimum 
cr(nt) — 0.03; however, constraints at this level are of interest for testing modifications to 
the the standard consistency relation in multifield inflation [101]. Constraints on Ug and 
Og do not improve significantly between the pre- and post-2020 timeframes, with errors de¬ 
creasing by a factor of approximately two in both cases. The minimum error on the energy 
density in curvature is (T(flk) = 1.5 x 10“^ (for Stage-IV x COrE-|-), five times stronger 
than current constraints but falling short of the 10““^ level interesting for testing the frame¬ 
work of eternal inflation when combined to other observations [102-107]. When considering 
the more challenging case of htting for the spatial variation in foreground spectral indices 
(in which case C-BASS or QUIJOTE-CMB are added as synchrotron monitors), A < 9.8 
in all cases apart from Stage-IV x Planck, with residuals remaining sub-dominant to the 
fiducial tensor signal. Consequently, a 3 - 7 -cj measurement of r = 10“^ is still attainable 
(after delensing) for certain combinations of experiments, even under the most challenging 
foreground scenario considered here. 

In the neutrino sector, the combination of Stage-IV with any space mission leads 
to constraints on the sum of the neutrino masses down to a{Mu) ~ 30 — 40 meV, and 
constraints on the effective number of relativistic species of a{Ng^) < 0.046, even under the 
A-expansion approach. Such measurements will begin to discriminate between the normal 
and inverted neutrino-mass hierarchies, and could confirm or spectacularly contradict our 
understanding of the thermal physics and neutrino makeup of the early universe [85]. The 
constraints on even the simplest model for the dark energy equation of state are relatively 
weak, with a{w) > 0.09. While this is considerably larger than the 1% constraints expected 
for Euclid [108] and LSST [67], adding a precision measurement of the Hubble parameter 
has the potential to significantly improve this constraint. Allowing a redshift-dependent 
dark energy equation of state yields errors on wo and Wa of ~ 0.10 and ~ 0.14 at best, 
though these numbers rely on delensing using measurements of large-scale structure: with 
access to only CMB and CIB data, these errors expand to ~ 0.25 and ~ 0.50, respectively. 

Overall, there is a clear synergy between the ground-based Stage-IV instruments 
(probing small angular scales, with deep integration at the CMB frequencies) and the 
space missions (probing large angular scales with broad frequency coverage). 

4 Web Interface 

In order to aid the community’s planning for future experimental designs, we make the 
toolkit used to obtain these results publicly available. The web interface is accessible at 
turkey.lbl.gov. The user is asked to fill in the following instrumental specifications: 

• central frequencies in GHz, 

• sensitivities in /rK-arcmin per frequency channel. 
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FWHM in arcmin per frequency channel, 

/sky) 

fractional bandpass per frequency channel, 


• ^min Q-nu tmax; 

• information channels among T, B and d - this is used to construct the covariance 
matrix, Eq. (2.29). 

In addition, the user can choose the foreground properties and analysis techniques, in 
particular: 

• dust temperature Tj and spectral index /3d, synchrotron spectral index Ps, 

• type of foreground separation among options described in Sect. 2.1, 

• type of delensing among the different options described in Sect. 2.2. 

• prior on /3d and /3s, Eq. (2.19) 

• the combination of the chosen specifications with Planck, following the procedure 
described in paragraph 2.5. 

Finally, the user has to choose the set of cosmological parameters for forecasts to be pro¬ 
duced. When these settings are chosen, the interface produces the following outputs: 

• the boosted noise in the reconstructed CMB map, ucmb, Eq- (2.16), 

• the level of the foreground residual power spectrum reff, Eq. (2.11), 

• the delensing factor a, Eq. (2.27), 

• the forecasted marginalized constraints on the chosen set of cosmological parameters. 

Wherever foregrounds are included in the input components, the constraints on cosmologi¬ 
cal parameters are also marginalized over the amplitude and tilt of the foreground residuals, 
cf. Sect. 2.3. 

5 Conclusions 

We have introduced a powerful and versatile framework for self-consistently evaluating the 
performance of next-generation CMB polarization experiments in the presence of astrophys- 
ical foregrounds and lensing. Our method, described in detail in Sect. 2, uses a parametric 
maximum-likelihood approach for component separation, CMB four-point estimators and 
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Figure 13. Top panel: Figure summarizing the performance of typical pre-2020 balloon and 
ground-based instruments, as well as their combination. The delensed B modes are obtained using 
CMB X CIB delensing. Bottom panel: Same as top panel but for post-2020 ground and space 
projects. Delensing is this time derived using the CMB x CMB estimator. 


cross-correlations with the cosmic infrared background and large-scale structure for delens¬ 
ing and the Fisher matrix formalism for parameter forecasts. The results of applying this 
framework to a suite of future ground-based, balloon and space-borne CMB polarization 
instruments were presented in Sect. 3. 
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Figure 13 provides a summary of the capabilities of this framework, demonstrating 
the ability of different experimental combinations to remove polarized foregrounds and 
lensing to reveal the primordial CMB. The frequency coverage and sensitivity of pre-2020 
experiments (top row) means they are noise-dominated in the search for B modes once 
foregrounds have been removed. Consequently, delensing does not improve constraints 
on primordial tensors. The post-2020 projects considered (bottom row) have much lower 
post-component-separation noise, even in the case of complex foreground emission, and 
their constraints on tensors are hence (de)lensing-limited. However, they will potentially 
be biased on large angular scales by foreground residuals. It is therefore important to note 
the power of the Stage-IV x space mission combination, which blends small angular scales 
and broad frequency coverage. 

Component separation — The component separation technique employed in this 
framework, based on a parametric maximum-likelihood method [47], consists of two steps. 
The first step is the estimation of the mixing matrix governing the frequency behaviour of 
the signal components, or more precisely the spectral indices of the polarized foregrounds; 
the second step is the linear combination of all available frequency maps to recover the 
cleaned CMB. Imperfect estimation of the mixing matrix leads to residual foregrounds 
contaminating the cleaned CMB maps, and inverting the mixing matrix to extract the 
CMB boosts the noise in the resulting maps. These effects are passed on to later stages 
of the pipeline, and the amplitude and angular dependence of the foreground residuals are 
marginalized over in all cosmological parameter constraints reported. 

Under the assumption that the foreground spectral indices are constant in independent 
patches of angular size ~ 15 deg, we show that Planck data are required for pre-2020 
projects to clean polarized dust and synchrotron, with the post-cleaning noise power spectra 
of the best experimental combinations approximately twice that expected from a naive 
combination of input frequency channels. Typical foreground residual amplitudes for these 
combinations correspond to an effective tensor-to-scalar ratio of reg < 5 x 10“^. As a 
result, the search for primordial B modes in this timeframe will be noise-dominated. Post- 
2020 instruments are more robust to the foreground composition, with noise degradations 
of A ~ 1.2 — 2.5 after component separation and foreground residual levels of rgg < 10“^. 
These residuals are comparable to the instrumental noise at the largest angular scales, and 
must therefore be carefully treated; the noise itself is sub-dominant to the lensing, and 
thus we expect delensing to be a powerful tool in the post-2020 timeframe. 

In an alternative approach to handling spatially varying spectral indices, we expand 
the mixing matrix to first order in the indices to allow them to vary from pixel to pixel. This 
approach is much more challenging, with noise degradations of A ~ 1.4 — 25.0 for even the 
post-2020 experiments, and requires dedicated synchrotron monitors such as C-BASS and 
QUIJOTE-CMB in addition to frequency bands above ~ 300 GHz in order to maximize 
the scientific output of the experiments considered. As emphasized in the A-expansion 
approach, C-BASS and QUIJOTE-CMB prove essential in disentangling the components 
in the case of an expansion of the mixing matrix, and the high-frequency coverage of PIXIE 
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makes it particularly robust in this more challenging scenario. This result, arising from the 
addition of new degrees of freedom to the mixing matrix, may hold in the cases of more 
complex foreground emissions (i.e., more components or more complex parametrizations). 

Exploration of the possible effects of mis-calibration or frequency band mismatch — 
potentially important systematic effects for component separation [47, 49] — is left for 
further work. 

Delensing — Our delensing method allows estimation of the lensing potential using 
the CMB polarization four-point function or the correlation between the CMB and two “ex¬ 
ternal” tracers of the lensing kernel: the GIB and large-scale structure. Delensing through 
cross-correlation, which requires only high-precision E mode measurements, outperforms 
“internal” CMB delensing for post-component-separation noise levels > 3 ^uK-arcmin, and 
represents the most promising avenue for pre-2020 experiments. We find that delensing 
using existing GIB data from Planck’’s 545 GHz channel is just as effective as using the 
highly futuristic large-scale structure measurements considered here, with both methods 
able to remove a maximum of 60% of the lensing using cosmic-variance-limited E-mode 
measurements (post-component-separation noise levels of < 10 /xK-arcmin). 

For the post-2020 experiments whose post-component-separation noise levels and res¬ 
olution allow precise measurements of the lensing B mode, the iterability of pure-CMB de¬ 
lensing means it rapidly outperforms the cross-correlation methods. In the limit of noiseless 
polarization measurements perfect delensing should be possible [29]; in the cases consid¬ 
ered here, up to around 80% of the lensing B modes can be removed by the combination 
of the ground-based Stage-IV and any satellite, assuming polarized dust and synchrotron 
have previously been cleaned from the map. Note that we do not explore in this work 
the possible impact of instrumental systematics which could limit lensing reconstruction, 
such as pointing mis-calibration, beam mis-characterization, optical systematics, etc. (see, 
e.g.. Ref. [109]); the impact of such instrumental systematics will grow with ever-increasing 
instrumental sensitivity. 

Parameter constraints — The Fisher estimate combines the two previous algo¬ 
rithms and produces constraints on any cosmological parameter. Our results are focused 
on r, TT-t, Us, Os, Mj,, w, and We show that in the case of GMB, synchrotron 
and dust, and after delensing and marginalization over foreground residuals, the best pre- 
2020 instruments in combination with Planck can reach (j{r) ~ 3 x 10“^, cj(nt) ~ 0.2, 
o'(ns) ~ 2.2 X 10“^, cr(as) ~ 3 x 10“^, a{M^) ~ 55 meV, a{w) ~ 0.16, ct(wo) ~ 0.36, 
o'iwa) ~ 0.71, (T(Neff) ~ 0.05 — 0.06 and (7(0^) ~ 2.5 x 10“^ when delensing using 
the GMB X GIB method. Post-2020 instruments, in particular the combination of the 
ground-based Stage-IV and a space mission, could reach constraints a{r) ~ 1.3 x lO”'^, 
<7(nt) ~ 0.03, cr(ns) ~ 1.8 x 10“^, (T(as) ~ 1.7 x 10“^, a{Mu) ~ 31 meV, a{w) ~ 0.09, 
(T(t(;o) ~ a{wa) ~ 0.50, cj(Vefj) ~ 0.024 and (7(11^) ~ 1.5x10“^. These constraints 
are robust under both approaches to spatially-varying foreground spectral indices. The 
combination of GMB measurements with Euclid, DESI, LSST and supernovae data will 
significantly improve the constraints on the dark sector. For example, constraints on w 
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derived from CMB observatories seem to be around 10%, compared with the target con¬ 
straints from Euclid or LSST {cr{w) ~ 1%). Finally, none of the CMB experiments studied 
is able to verify the inflationary consistency condition nt = —r/8 after disentangling dust, 
synchrotron and CMB, even assuming a level r = 0.1 of primordial B modes (confirming 
the results of Ref. [41], which are derived in a more optimistic setting). Again, external 
priors like H{z) given by e.g., DESI could significantly improve all the quoted marginalized 
constraints. 

Finally, we have made our forecasting toolkit publicly available through a web in¬ 
terface at turkey. Ibl. gov. We hope that this tool and its extensions will serve the 
community by enabling the optimal design of future CMB observatories for a broad range 
of scientific goals, thus unleashing the full potential of the polarized CMB sky to constrain 
fundamental physics. 
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A Instrument specifications 

Tables 4 to 6 give the specifications of each instrument considered in this work. The 
maximum multipole tabulated for each instrument is indicative: as in Fig. 5, it is taken to 
be the minimum of 4000 and the multipole at which the beam-deconvolved noise reaches 
10^where it dominates over even the temperature power spectrum. Fig. 15 shows 
polarized sensitivities of all post-2020 instruments as a function of frequency. 
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Table 4. Planck, C-BASS, QUIJOTE-CMB 


Planck specifications, Ref. [23 


frequencies [GHz] 

fractional bandpass [%] 

sensitivities [^K-arcinin] 

/sky [%] 

FWHM [arcinin] 

P • 

^min 

P 

'-max 

30.0 


300.0 


33.2 



44.0 


300.0 


28.1 



70.0 


192.0 


13.1 



100.0 

30.0 

52.0 

50.0 

9.7 

3 

2800 

143.0 


44.0 


7.3 



217.0 


64.0 


5.0 



353.0 


276.0 


4.9 



C-BASS specifications, http://www.astro.caltech.edu/cbass 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [/iK-arcmin] 

/sky [%] 

FWHM [arcinin] 

P • 

^min 

P 

'-max 

5.0 

20.0 

4500.0 

80.0 

45.0 

20 

500 

QUIJOTE-CMB specifications, Ref. [90] 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [/iK-arcmin] 

/sky [%] 

FWHM [arcinin] 

P • 

^min 

P 

'-max 

11.0 

18.0 

840.0 


55.2 



13.0 

15.0 

840.0 


55.2 



17.0 

12.0 

840.0 

44.0 

36.0 

20 

500 

19.0 

11.0 

840.0 

36.0 

30.0 

27.0 

66.6 


22.2 



42.0 

24.0 

50.4 


16.8 





Figure 14. Polarized sensitivities for PIXIE. Bandpasses are 15 GHz wide across all frequency 
channels. 
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Table 5. Pre-2020 instruments 


Advanced ACTPol specifications, Ref. [78] 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [//K-arcmin] 

Aky [%] 

FWHM [arcinin] 

^min 

^max 

90.0 


11.0 


2.2 



150.0 

30.0 

9.8 

50.0 

1.3 

20 

4000 

230.0 


35.4 


0.9 



BICEP3 + Keck specifications 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [//K-arcmin] 

/sky [%] 

FWHM [arcinin] 

^min 

^max 

95.0 

30.0 

1.7 

1.0 

25.0 

20 

1300 

150.0 

3.4 

30.0 

CLASS specifications, Ref. [80 


frequencies [GHz] 

fractional bandpass [%] 

sensitivities [//K-arcmin] 

/sky [%] 

FWHM [arcinin] 

^min 

^max 

38.0 

93.0 

30.0 

39.0 

10.0 

70.0 

90.0 

40.0 

20 

1100 

148.0 

15.0 

24.0 

217.0 


43.0 


18.0 



EBEXlOK specifications, proposal to NASA in 2015 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [//K-arcmin] 

/sky [%] 

FWHM [arcinin] 

^min 

^max 

150.0 

220.0 

30.0 

5.5 

11.0 

2.5 

6.6 

4.7 

20 

4000 

280.0 

25.4 

3.9 

350.0 


53.0 


3.3 



PIPER specifications, Ref. [83 


frequencies [GHz] 

fractional bandpass [%] 

sensitivities [^K-arcmin] 

/sky [%] 

FWHM [arcinin] 

^min 

^max 

200.0 

30.0 

31.4 


21.0 



270.0 

30.0 

45.9 

85.0 

21.0 

20 

1000 

350.0 

16.0 

162.0 

21.0 

600.0 

10.0 

2659.2 


21.0 



Simons Array specifications, Ref. [82] 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [/./K-arcinin] 

/sky [%] 

FWHM [arcmin] 

f ■ 

^min 

f 

^max 

90.0 


14.4 


5.2 



150.0 

30.0 

11.8 

65.0 

3.5 

20 

4000 

220.0 


40.3 


2.7 



SPIDER specifications, Ref. [84] 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [/./K-arcinin] 

/sky [%] 

FWHM [arcinin] 

f ■ 

f 

^max 

90.0 

24.0 

21.2 

8.0 

45.0 

20 

800 

150.0 

24.0 

17.7 

30.0 

SPT-3G specifications, Ref. [81] 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [//K-arcinin] 

/sky [^1 

FWHM [arcmin] 

f ■ 

P 

^max 

95.0 

27.0 

7.0 


1.6 



148.0 

26.0 

4.5 

6.0 

1.1 

20 

4000 

223.0 

23.0 

7.5 


1.0 
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Table 6. Post-2020 instruments 



GOrE+ specifications, http://conservancy 

urDn.edu/handle/11299/169642 


frequencies [GHz] 

fractional bandpass [%] 

sensitivities [^tK-arcinin] 

/at [%] 

FWHM [arcmin] 

^inin 


60.0 


16.0 


14.0 



70.0 


14.9 


12.0 



80.0 


12.9 


10.5 



90.0 


9.2 


9.3 



100.0 


8.5 


8.4 



115.0 


7.0 


7.3 



130.0 


5.9 


6.5 



145.0 


5.0 


5.8 




30.0 


70.0 


2 


160.0 


5.4 


5.3 



175.0 


5.3 


4.8 


4000 

195.0 


5.3 


4.3 



220.0 


8.1 


3.8 



255.0 


12.6 


3.3 



295.0 


27.4 


2.9 



340.0 


43.7 


2.5 



390.0 


77.8 


2.2 



450.0 


164.8 


1.9 



520.0 


418.2 


1.6 



600.0 

_ 

1272.4 


1.4 



LiteBIRD-ext specifications http: //ltdl6. grenoble. cnrs. fr/IMG/UserFiles/Images/09_TMatsumura_20150720_LTD_vl8. pdf 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [/rK-arcmin] 

/sify [%] 

FWHM [arcmin] 

^miii 


40.0 


42.5 


108 



50.0 


26.0 


86 



60.0 


20.0 


72 



68.4 


15.5 


63 



78.0 


12.5 


55 



88.5 


10.0 


49 



100.0 


12.0 


43 



118.9 

30.0 

9.5 

70.0 

36 

2 

1350 

140.0 


7.5 


31 



166.0 


7.0 


26 



195.0 


5.0 


22 



234.9 


6.5 


18 



280.0 


10.0 


37 



337.4 


10.0 


31 



402.1 


19.0 


26 



Stage-IV specifications, derived so that the noise after component separation, (tcmb- is ~ 1 pK-arcmin, Refs. [70, 85] 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [/rK-arcmin] 

/sisy [%] 

FWHM [arcmin] 


^max 

40.0 


3.0 


11.0 



90.0 


1.5 


5.0 



150.0 

30.0 

1.5 

50.0 

3.0 

20 

4000 

220.0 


5.0 


2.0 



280.0 


9.0 


1.5 



PIXIE specifications, Ref. [86] 

frequencies [GHz] 

fractional bandpass [%] 

sensitivities [/^K-arcmin] 

/.ky [%I 

FWHM [arcmin] 

^inin 

^max 

see Fig. 14 

70.0 

96.0 

2 

500 
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sensitivity sensitivity sensitivity sensitivity 
[/iK-arcmin] [/iK-arcmin] [/iK-arcmin] [/iK-arcmin] 
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Figure 15. Comparison of polarized sensitivities for post-2020 instruments between 30 and 800 
GHz. 




































































































B Error bars on cosmological parameters for each project 

We detail in this section the error bars on several relevant cosmological parameters as esti¬ 
mated by our algorithm, for the cross-correlation of all pre-2020 instruments with Planck, 
as well as for the post-2020 observatories. In addition, we summarize the delensing factor, 
the effective level of foreground residuals and the noise degradation introduced in Sect. 2. 
Each constraint is derived assuming a particular cosmology, as detailed in Table 3. 

In the cases considering combinations of datasets, A and reff are the values obtained 
over the overlapping sky, i.e., when all the frequency maps and sensitivities from the two 
instruments are combined. 

B.l Pre-2020 instruments, cross-correlated with Planck 


BICEP3 X Keck Array x Planck 

Delensing option —> 

no 

CMB X 

CMB 

X 

CMB 

X 

comp. sep. option 

delensing 

CMB 

CIB 


LSS 




a = 

1.0 

a = 

0.66 

a = 0.43 

a = 0.43 


A = 1.0 

o-(r) = 2.5 X 

a{AU) = 1.7 X 10® 

iT{r) = 2.1 X 10“® 

tr(M„) = 1.7 X 10® 

o-(r) = 1.6 X 10“® 

tr{M„) = 69 

tr(r) = 1.6 X 10“® 

(t(A4) = 50 



cr(nt) = 0.10 

cr(Weff) = 0.15 

(7{nt) = 0.098 

fr(Aeff) = 0.15 

cr(nt) = 0.094 

0-{Aeff) = 0.10 

= 0.094 

o-Aeff) = 0.069 

cln oil y 


a{n,) = 3.0 X 10“® 

a{w) = 0.41 

(7(ns) = 3.0 X 10“^ 

cr{w) = 0.41 

cr(ns) = 3.0 X 10-® 

= 0.18 

tr(ne) = 2.9 X 10-® 

(T(itf) = 0.13 


j-eif = o.n 

o-(as) = 4.8 X 10“® 

cr(ii;o) = 1.3 

cr{as) = 4.8 X 10“^ 

ir(mo) = 1.3 

cr(as) = 4.4 x 10“^ 

cr{wQ) = 0.45 

o-(as) = 2.6 X 10“® 

(t(ico) = 0.14 



o-(nt) = 7.0 X 10-® 

a[wa) = 2.6 

D-iflk) = 7.0 X 10-® 

a{wa) = 2.6 

= 2.9 X 10-® 

= 0.87 

o-(f!t) = 2.6 X 10-® 

cr{wa) = 0.23 



a = 

1.0 

a = 

0.97 

a = 0.53 

a = 0.54 


A = 53 

o-(r) = 0.013 

tr(M„) = 3.0 X 10® 

iT{r) = 0.013 

tr(M„) = 3.0 X 10® 

a{r) = 0.013 

= 74 

o-(r) = 0.013 

= 52 



cr(nt) = 0.23 

cr(Weff) = 0.16 

(7{nt) = 0.23 

fr(Aeff) = 0.16 

cr(nt) = 0.23 

0-{Aeff) = 0.11 

cr(nt) = 0.22 

o-Aeff) = 0.073 

sylic US 


cr{n,) = 3.3 x 10“® 

fr(w) = 0.61 

a(nj) = 3.3 x 10“® 

a{w) = 0.61 

cr(ns) = 3.2 X 10-® 

(7{it;) = 0.19 

ir(iie) = 3.2 X 10-® 

= 0.13 


reff = 2.7 X 10“® 

o-(as) = 5.3 X 10^4 

cr(tco) = 1.4 

cr{as) = 5.3 x 10“^ 

tr(mo) = 1.4 

cr(as) = 4.7 x 10“^ 

= 0.45 

tr(as) = 2.7 x 10“® 

cr('ii;o) = 0.14 



0-(S4) = 0.014 

a{wa) = 3.1 

(7{nk) = 0.014 

a{Wa) = 3.1 

= 3.1 X 10-® 

a{wa) = 0.89 

tr(rik) = 2.8 X 10-® 

cr[wa) = 0.24 


EBEXlOK X Planck 

Delensing option —)• 

4, comp. sep. option 

no 

delensing 

CMB X 

CMB 

CMB X 

CIB 

CMB X 

LSS 

cmb-oiily 

A = 1.0 

reff = 0.0 

a = 

a{r) = 2.7 X 10-® 

CT(rit) = 0.10 

ct(iIj) = 2.9 X 10-® 
CT(as) = 4.6 X 10-® 

cr(nk) = 5.1 X 10-® 

1.0 

a{AU) = 1.2 X 10^ 
CT(Aeff) = 0.13 
ct(u;) = 0.33 

CT(n;o) = 1.2 

a{Wa) = 2.4 

a = 

cr(r) = 2.5 X 10-® 
CT(ni) = 0.099 
ffK) = 2.9 X 10-® 
ff(Qs) = 4.6 X 10-® 
ffifik) = 5.1 X 10-® 

0.61 

(t{AU) = 1.2 X 10^ 
o-(Ae£f) = 0.13 
crfic) = 0.33 
cr(ico) = 1.2 
a{wa) = 2.4 

a = 0 

o-(r) = 2.3 X 10-® 
o-(ni) = 0.097 
o-(n,) = 2.9 X 10-® 
o-(Qs) = 4.3 X 10-® 
o-(Slk) = 2.8 X 10-® 

41 

cr(M^) = 66 
o-(Aeff) = 0.099 
cr(w) = 0.17 
tT(n;o) = 0.45 

o-(wa) = 0.86 

a = 0 

CT(r) = 2.3 X 10-® 
o-(ni) = 0.097 
a(n/) = 2.8 x 10-® 
<t(Qs) = 2.6 X 10-® 
ffifik) = 2.5 X 10-® 

39 

a(A/^) = 48 
a(JVeB) = 0.067 
a(w) = 0.13 
(t(ico) = 0.14 
fr(w„) = 0.23 

sync-hdust 

A = 2.2 

reff = 8.8 X 10-^ 

a = 

CT(r) = 6.6 X 10-® 
a(nt) = 0.21 
cr(ns) = 3.1 X 10-® 
cr(ae) = 5.0 X 10-® 
crifik) = 6.9 X 10-® 

1.0 

a{M^) = 1.7 X 10^ 
a(Aeff) = 0.15 
cr{w) = 0.41 
cr{wQ) = 1.3 
a{wa) = 2.7 

a = 

ff(r) = 6.1 X 10-® 
cr(nt) = 0.21 
cr(ne) = 3.1 X 10-® 
cr(Qs) = 4.9 X 10-® 
crifik) = 6.9 X 10-® 

0.73 

cr(A4) = 1.7 X 10^ 
'r(77eff) = 0.15 
a{w) = 0.41 
o-(n;o) = 1.3 
cr(tf;Q) = 2.7 

a = 0 

a(r) = 5.5 X 10-® 
a-(nt) = 0.21 
tr(ns) = 3.1 X 10-® 
tr(as) = 4.5 X 10-® 
o-(Slk) = 3.0 X 10-® 

41 

d-(M^) = 72 

cr(Aeff) = 0.11 

o-(w) = 0.18 
cr(n;o) = 0.45 
cr(wa) = 0.88 

a = 0 

a(r) = 5.5 x 10“^ 

(T(nt) = 0.20 

cr(ns) = 3.0 X 10"^ 
cr(as) = 2.7 X 10"^ 
cr(Dk) = 2.7 X 10"^ 

40 

cr(A4) = 51 
a(JV^ff) = 0.070 
a(w) = 0.13 
a(wo) = 0.14 
a'(Wa) = 0.24 


- 42 - 




























































AdvACTPol X Planck 

Delensing option 

1 comp. sep. option 

no 

delensing 

CMB X 

CMB 

CMB X 

CIB 

CMB X 

LSS 

cmb-only 

A - 1.0 

TeS - 0.0 

Q = 1 

tT(r) = 1.3 X 10“^ 
cT(jit) — 0.087 
alris) = 2.2 x 10“=‘ 
tT(as) = 2.6 X 10“^ 
t7(fik) = 2.2 X 10“^ 

.0 

a{My) — 53 
C7(Aeff) - 0.044 
a{w) = 0.14 
(7{wq) = 0.51 
a{wa) = 1.0 

a = 0. 

(j{r) = 1.2 X 10“^ 
fT(jit) — 0.084 
(T(raa) = 2.2 X 10“^ 
(T(aa) = 2.6 X 10“** 
(T(rik) = 2.2 X 10“^ 

66 

a[My) — 53 
a{N^s) - 0.044 
a{w) = 0.14 
cT(it;o) = 0.51 
a{wa) = 1.0 

a = 0. 

(T(r) — 1.1 X 10“^ 

fT(jit) — 0.082 

cr(ns) - 2.2 X 10“3 
(T(as) - 2.6 X 10“^ 
(T(f^k) - 2.1 X 10“^ 

41 

(7(Mi,) — 49 
a(N,ff) - 0.042 
(j(w) = 0.14 
(t(i[;o) — 0.36 
— 0.67 

a = 0. 

(T(r) = 1.1 X 10“® 

fT(nt) — 0.081 

r7(nj = 2.2 x 10“^ 
(T(Qa) = 2.1 X 10“^ 
(T(rik) = 2.0 X 10“^ 

39 

^ 39 

a(N^ff) - 0.034 
(7(w) — 0.11 
(t(ujo) — 0.12 
— 0.18 

sync+dust 

A - 3.6 

^’eff — 3.8 X 10“^ 

Q — 1 

rj{r) = 3.8 x 10“^ 

CT(jlt) — 0.16 

tj{ns) = 2.3 X lO”** 
tj{as) = 3.1 X 10“^ 
tT(fik) = 3.3 X 10“=* 

.0 

a{M^) - 78 

C7(Aeff) - 0.060 

a[w) — 0.19 
C7(lI!o) — 0.69 
a{wa) = 1.4 

a = 0. 

(j{r) = 3.7 X 10“^ 

(T{nt) — 0.16 

<T(raa) = 2.3 X 10“^ 
(T(aa) = 3.1 X lO”** 
(T(fik) = 3.3 X 10“^ 

85 

a{M^) - 78 

C7(AefF) - 0.060 

a{w) — 0.19 
ct(ii!o) — 0.69 
a{wa) = 1.4 

a = 0. 

cr(r) = 3.5 X 10“^ 

(T(jit) — 0.16 

- 2.3 X 10“3 
(T(as) - 3.0 X 10“^ 
(T(f)k) = 2.4 X 10“^ 

43 

o-(Mi,) — 55 
a(N,ff) - 0.054 

(j(w) = 0.16 

(t(ii!o) — 0.40 
(7(Wa) = 0.74 

a = 0. 

cr(r) = 3.5 X 10“^ 
(T(rtt) — 0.15 
fT(7is) - 2.3 X 10“^ 
(t(qs) - 2.3 X 10“^ 
(7(fik) = 2.2 X 10"^ 

41 

a(My) - 42 
^(N,ff) - 0.042 

(7(w) = 0.12 

(t(iuo) — 0.13 

(7(Wa) = 0.20 


PIPER X Planck 

Delensing option — 

no 

CMB X 

CMB 

X 

CMB 

X 

1 comp. sep. option 

delensing 

CMB 

CIB 


LSS 




a = 

1.0 

a = 

0.97 

1 a = 0.50 1 

1 a = 0.50 1 


A = 1.0 

a{r) = 3.8 X 10-^ 

cr(M^) = 2.3 X 10^ 

a{r) = 3.8 x 10“^ 

a{M^) = 2.3 X 10^ 

o-(r) = 3.7 X 10-3 

cr(M^) = 59 

<T(r) = 3.7 X 10-3 

a{M^) = 40 

cmb-only 


(r(nt) = 0.12 
<T(rts) = 2.4 X 10-® 

o-(Aeff) = 0.13 
o-(tt)) = 0.51 

(r{nt) = 0.12 

a{nj) = 2.4 x 10“^ 

<T(Ale£l) = 0.13 
<T(tt;) = 0.51 

cr(nt) = 0.12 
o-(iis) = 2.4 X 10-3 

o-(Aeff) = 0.088 

cr(ic) = 0.15 

= 0.12 

ff(iis) = 2.4 X 10-3 

o-(Aeff) = 0.058 
a{w) = 0.11 


reff = 0.0 

CT(aa) = 4.1 X 10“^ 

o-(tt)o) = 1.2 

a(a^) = 4.1 X 10“^ 

(t(ico) = 1-2 

o-ftta) = 3.7 X 10-3 

cr(ii)o) = 0.36 

(t(q:s) = 2.2 X 10-3 

o-(ii)o) = 0.11 



ff(!lk) = 0.010 

O'(tCa) = 2.6 

o-{nk) = 0.010 

a{wa) = 2.6 

o-(r!k) = 2.5 X 10-3 

<7(tCa) = 0.69 

(T(f!k) = 2.2 X 10-3 

cr{Wa) = 0.19 



a = 

1.0 

a = 

0.98 

1 a = 0.51 1 

1 a = 0.51 1 


A = 1.2 

a{r) = 8.3 x 10“^ 

o-(M^) = 2.6 X 10^ 

cr{r) = 8.3 x 10"^ 

a{M^) = 2.6 X 10^ 

o-(r) = 8.1 X 10-3 

o-(M^) = 68 

<T(r) = 8.1 X 10-3 

a{M^) = 44 

sync-|-diist 


a{nt) = 0.20 
ffK) = 2.8 X 10-® 

O'(A^eff) = 0.15 

D-(ic) = 0.59 

o‘{nt) = 0.20 
cr{ns) = 2.8 X 10"^ 

<T(A^eff) = 0.15 
<T(tt;) = 0.59 

o-(nt) = 0.20 
o-(iis) = 2.8 X 10-3 

o-(Aeff) = 0.093 

o-(tt)) = 0.16 

a(nt) = 0.20 
<t(iIs) = 2.8 X 10-3 

o-(Aeff) = 0.060 
a{w) = 0.11 


reff = 4.1 X 10“^ 

fffas) = 4.8 X 10“^ 

(7(tt)o) = 1.4 

cr{as) = 4.8 X 10“^ 

a{wo) = 1.4 

o-(qs) = 4.1 X 10-3 

a{w(j) = 0.36 

(t(q:s) = 2.3 X 10-3 

cr(tt)o) = 0.12 



<T(!2k) = 0.012 

(7{'Wa) = 2.9 

cr{Dk) = 0.012 

cr{Wa) = 2.9 

irislk) = 2.9 X 10-3 

O'(tCa) = 0.71 

ffislk) = 2.5 X 10-3 

cr{Wa) = 0.19 


1 SPIDER X Planck \ 

Delensing option —> 

1 comp. sep. option 

no 

delensing 

CMB X 

CMB 

CMB X 

CIB 

CMB X 

LSS 

cmb-only 

A = 1.0 

I'eff = 0.0 

a = 

(7{r) = 3.1 X 10“^ 
fr{nt) = 0.11 
cT{ns) = 3.0 X 10-3 
(7{as) = 4.9 X 10-3 
o-{nk) =0.011 

1.0 

o-(M^) = 2.6 X 10^ 
fr(Weff) = 0.15 
o-(tt)) = 0.52 
o-(wjo) = 1.3 
a{wa) = 2.8 

a = 

<T(r) = 3.1 X 10-3 
= 0.11 

tT(n,) = 3.0 X 10-3 
(t(qs) = 4.9 X 10-3 
<T(r!t) = 0.011 

0.97 

a{M^) = 2.6 X 10^ 
a{Nes) = 0.15 
a{w) = 0.52 
(t(i(Jo) = 1.3 
a{wa) = 2.8 

a = 0.^ 

(7{r) = 3.1 X 10-3 

cT(nt) = 0.11 
(7(ns) = 2.9 X 10-3 
(7(as) = 4.4 X 10-3 
c7(£2k) = 2.9 X 10-3 

19 

o-{M„) = 69 
= 0.10 
o-{m) = 0.18 
fT{w;o) = 0.45 
tj{w^) = 0.87 

Q = 0. 

cr(r) = 3.1 X 10-3 

o-{nt) = 0.11 

tr(ns) = 2.9 x 10-3 
(T(Qa) = 2.6 X 10-3 
tr(nk) = 2.6 X 10-3 

50 

= 50 

(T(Aeff) = 0.069 
a{w) = 0.13 
(T(t(;o) = 0.14 
a{wa) = 0.23 

sync-|-dust 

A = 2.2 

reff = 4.8 X 10-^ 

a = 

a(r) = 0.012 
fT{nt) = 0.22 

(7{ns) = 3.2 X 10-3 
(7{as) = 5.3 X 10-3 
cT{r2k) = 0.014 

1.0 

o-(M^) = 3.1 X 10^ 
f’‘(A^eff) = 0.16 
(t{w) = 0.63 
a-(wjo) = 1.4 
a{wa) = 3.1 

a = 

a{r) = 0.012 
(T(nt) = 0.22 

a{ns) = 3.2 x 10-3 
(t(q:s) = 5.3 X 10-3 
(T(f2k) = 0.014 

0.99 

a{M^) = 3.1 X 10^ 

a(Aeff) = 0.16 

= 0.63 
<t(h;o) = 1.4 
<T(u)a) = 3.1 

a = 0.' 

(7[r) = 0.012 
cr(nt) = 0.22 
(7(ns) = 3.2 X 10-3 
(7(ag) = 4.7 X 10-3 
cT(r2k) = 3.1 X 10-3 

>6 

a{M^) = 74 

fT{Aeff)=0.11 

a{w) = 0.19 
<7{wq) = 0.45 
(T{Wa) = 0.89 

Q = 0. 

o-(r) = 0.012 

(T{nt) = 0.22 

tr(ns) = 3.2 x 10-3 
tr(aa) = 2.7 x 10-3 
trifik) = 2.8 X 10-3 

56 

a{A'U) = 52 
a{Nes) = 0.073 
(t(ju) = 0.13 
(t(j(;o) = 0.14 
a{wa) = 0.24 


Simons Array x Planck 

Delensing option 

no 


1 CMB 

X 

CMB 

X 

CMB 

X 

1 comp. sep. option 

delensing 

CMB 

CIB 


LSS 




Q = 1.0 

a = 0.76 

1 a = 0.41 1 

a = 0.40 


o 

II 

<1 

cT(r) — 1.6 X 10 3 

cT(Af^) = 56 

(T(r) — 1.5 X 10 3 

a{M^) = 56 

(T(r) — 1.4 X 10 3 

(7{M^) = 47 

(7{r) — 1.4 X 10 3 

a{M^) = 37 

cmb onl 


cT(jit) — 0.089 

o-(A^eff) - 0.055 

fT(jit) = 0.088 

c^(A'eff) — 0.055 

= 0.085 

fT(A^efF) - 0.049 

fT(nt) = 0.085 

(7{Nes) - 0.038 

cm -ony 


tjjraa) = 2.0 X 10-3 

a{w) = 0.14 

(T(raa) = 2.0 X 10-3 

a{w) = 0.14 

(T(raa) = 2.0 X 10-3 

a{w) = 0.14 

(T(na) = 2.0 X 10-3 

(7{w) — 0.11 


fell = 0.0 

tT(aa) = 2.7 X 10-3 

(j{wq) — 0.63 

(T(aa) = 2.7 X 10-3 

(j{wq) — 0.63 

(T(aa) = 2.7 X 10-3 

(T(r[;o) — 0.36 

(T(Qa) = 2.0 X 10-3 

(T(rtJo) — 0.11 



tj(fik) = 2.3 X 10-3 

= L3 

(T(fik) = 2.3 X 10-3 

cT(it;a) = L3 

(T(fik) = 2.0 X 10-3 

(7(it;a) = 0.66 

(T(rik) = 1.9 X 10-3 

(^{Wa) = 0.18 



1 Q=1.0 1 

1 a = 0.89 1 

1 a = 0.43 1 

1 a = 0.42 1 


A - 3.0 

c7(r) — 4.6 X 10 3 

a{M^) = 96 

(T(r) = 4.5 X 10 3 

a{M^) = 96 

(T(r) = 4.3 X 10-3 

= 55 

(T(r) = 4.3 X 10-3 

= 40 

sync-j-dust 


CT(jlt) — 0.16 

o-(A^eff) - 0.076 

fT(71t) = 0.16 

C7(AefF) = 0.076 

(T(jit) = 0.16 

(T(Aeff) = 0.064 

(T(nt) = 0.16 

(j{Nes) = 0.046 



tT(raa) = 2.2 X 10-3 

a{w) = 0.23 

(T(raa) = 2.2 X 10-3 

a{w) = 0.23 

(T(raa) = 2.2 X 10-3 

(7{w) = 0.16 

(7{nJ) = 2.2 X 10-3 

(7{w) — 0.11 


7-eff - 5.0 X 10-3 

t7(aa) = 3.3 X 10-3 

cT(it;o) = 0.84 

cr(aa) = 3.3 X 10-3 

(j{wq) — 0.84 

(T(aa) = 3.1 X 10-3 

(T(r[!o) — 0.38 

(j(a,) = 2.2 X 10-3 

(T(r[Jo) = 0.12 



t7(fik) =4.0 X 10-3 

II 

b 

(T(fik) = 4.0 X 10-3 

(^(w'a) = L7 

CT(fik) = 2.3 X 10-3 

= 0.71 

(T(rik) = 2.1 X 10-3 

'^{Wa) = 0.20 
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1 CLASS X Planck \ 

Delensing option —> 
i comp. sep. option 

no 

delensing 

CMB X 

CMB 

CMB X 

CIB 

CMB X 

LSS 

cmb-only 

A = 1.0 

refi = 0.0 

a = 

a{r) = 1.4 X 10-^ 
a{nt) = 0.092 
a{nj) = 2.4 x 10“® 
<T(as) = 4.2 X 10-® 
CT(nk) = 6.5 X 10“® 

1.0 

a{M^) = 1.6 X 10^ 
CT(Veff) = 0.13 
ct(w) = 0.35 

(7(«fo) = 1.2 

(7{Wa) = 2.5 

a = 1 

a(r) = 1.4 X 10-® 
cr(nt) = 0.092 
a\nt) = 2.4 x lO^® 
ala,) = 4.2 x 10“^ 
CT(!!k) = 6.5 X 10“® 

0.95 

cr(A/^) = 1.6 X 10^ 
a-(7Veff) = 0.13 
cr(t(f) = 0.35 
a-(«fo) = 1-2 
<7{Wa) = 2.5 

(1 = 0. 

a{r) = 1.2 X 10-® 
alrit) = 0.089 
a{n,) = 2.4 x 10“^ 
o-Ls) = 3.8 X 10-3 
i7(!ik) = 2.5 X 10-3 

48 

cr(M^) = 60 
a-(Aefl) = 0.092 
cr(w) = 0.16 

cr(wo) = 0.39 
a(wa) = 0.75 

(1 = 0. 

a(r) = 1.2 X 10-3 
alrtt) = 0.089 
aln,) = 2.4 x 10-3 
ala,) = 2.3 X 10-3 
cr(!!k) = 2.2 X 10-3 

48 

a(M^) = 42 
'^(A^eff) = 0.061 
a(w) = 0.11 
a(wo) = 0.12 

(7(tt;„) = 0.20 

sync+dust 

A = 2.1 

refi = 3.5 X 10“^ 

a. = 

a(r) = 3.0 X 10-^ 

= 0.15 

CT(ns) = 2.7 X 10-® 
<T(Qa) = 4.7 X 10“® 
<T(nk) = 0.011 

1.0 

a{M^) = 2.5 X 10^ 

= 0.15 
u{w) = 0.55 
a{wQ) = 1.4 
(7{Wa) = 2.9 

a = 1 

a(r) = 3.0 X 10-® 
a{n,) = 0.15 
ctK) = 2.7 X 10-^ 
ala,) = 4.7 X 10“^ 

iT(!2k) = 0.011 

0.98 

cr(i\4) = 2.5 X 10^ 
a{N^fi) = 0.15 
alw) = 0.55 
cr(tt;o) = 1.4 
a{Wa) = 2.9 

(1 = 0. 

a{r) = 2.8 x 10-3 
a{nt) = 0.15 
aln,) = 2.7 X 10-3 
ala,) = 4.1 X 10-3 
i7(S2k) = 2.8 X 10-3 

52 

(t(M^) = 67 
O'(A^efl) = 0.10 
(t(w) = 0.17 
cr(n;o) = 0.40 
a(wa) = 0.77 

(1 = 0. 

a(r) = 2.8 X 10-3 
a(nt) = 0.15 
ctK) = 2.6 X 10-3 
a(as) = 2.4 x 10-3 
<T(!2k) = 2.4 X 10-3 

53 

(7(i\4) = 45 

= 0.066 
a(w) = 0.12 
a{wQ) = 0.12 
a{Wa) = 0.21 


SPT-3G X Planck 

Delensing option —> 
i comp. sep. option 

no 

delensing 

CMB X 

CMB 

CMB X 

CIB 

CMB X 

LSS 

cmb-only 

A = 1.0 

refi = 0.0 

a = 

a{r) = 1.7 X 10-3 
alnt) = 0.090 

17+) = 2.6 X 10-3 
ala,) = 3.7 x 10-3 
<T(nk) = 2.8 X 10-3 

1.0 

a{M^) = 66 
a(7Veff) = 0.073 
a{w) = 0.20 
a{wQ) = 0.76 
a{Wa) = 1.5 

a = 1 

a(r) = 1.1 X 10-3 
ct(tii) = 0.083 

<7+) = 2.6 X 10-3 
ala,) = 3.7 X 10-3 
i7(!2k) = 2.8 X 10-3 

0.40 

cr(M^) = 66 

= 0.073 
a{w) = 0.20 
a{wQ) = 0.76 
cr(tt;„) = 1.5 

(1 = 0. 

£7(r) = 1.1 X 10-3 

£7(nt) = 0.083 
a(n,) = 2.6 X 10-3 
£7(£i8) = 3.5 X 10-3 
£7(f!k) = 2.4 X 10-3 

40 

cr(M^) = 56 
o-(A^efi) = 0.065 

cr(w) = 0.16 

cr(n;o) = 0.41 
cr(Wa) = 0.79 

(1 = 0. 

£7(r) = 1.1 X 10-3 

alnt) = 0.083 
aln,) = 2.6 x 10-3 
ala,) = 2.4 x 10-3 
£r(f!k) = 2.2 X 10-3 

38 

a(M^) = 44 

= 0.050 
a(w) = 0.12 
a(wo) = 0.13 

a(Wa) = 0.21 

sync+dust 

A = 5.9 

refi = 5.7 X 10-^ 

a = 

a{r) = 4.8 x 10-3 
a{nt) = 0.19 
£7+) = 2.8 X 10-3 
ala,) = 4.1 X 10-3 

ffifik) = 4-2 X 10-3 

1.0 

a{M^) = 1.0 X 10^ 
<7(A^eff) = 0.097 
ct(w;) =0.27 

ct(w;o) = 1.0 
a{Wa) = 2.1 

a = 1 

<7(r) = 4.5 X 10-3 
a{nt) = 0.19 
£7+) = 2.8 X 10-3 
ala,) = 4.1 X 10-3 
£7(f!k) = 4.2 X 10-3 

0.70 

ct{Mu) = 1.0 X 10^ 
o-(-^eff) = 0.097 
a{w) = 0.27 
(T(-i(fo) = LO 

a{Wa) = 2.1 

(1 = 0. 

a(r) = 4.1 X 10-3 
(7(nt) = 0.19 
aln,) = 2.7 X 10-3 
£7(q:8) = 3.9 X 10-3 
£7(f!k) = 2.6 X 10-3 

41 

cr(M^) = 62 

o-(N^ff) = 0.080 

(T(n;) = 0.17 
cr(wo) = 0.44 
cr(Wa) = 0.83 

(1 = 0. 

£7(r) = 4.1 X 10-3 
a(nt) = 0.18 
£7+) = 2.7 X 10-3 
£7+) = 2.5 X 10-3 
£7(f!k) = 2.4 X 10-3 

39 

= 46 
= 0.058 
a(w) = 0.12 
a(wo) = 0.14 
(7(tt;„) = 0.23 
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B.2 Post-2020 instruments 


In addition to the Up approach used for the sync+dust cases shown in previous tables, 
the following results include the A-expansion approach after combining the considered 
instruments with either C-BASS or QUIJOTE-CMB. Independently of their specific sky 
coverage (80% for C-BASS and 44% for QUIJOTE-CMB), we assume these two obser¬ 
vatories help characterize synchrotron emission over the sky observed by any post-2020 
instrument. Moreover, QUIJOTE-CMB is only combined with space missions as it is a 
project observing the Northern hemisphere; as with all the ground-based Stage-III instru¬ 
ments we consider throughout this work, we assume that the Stage-IV project will likely be 
observing the Southern hemisphere. We add Planck when considering Stage-IV in isolation 
in order to include information about the largest angular scales. 


Stage-IV X Planck 

Delensing option — > jj no 

1 comp. sep. option delensing 

CMB X 

CMB 

CMB X 

CIB 

CMB X 

LSS 

cmb-only 

A - 1.0 

r-eff ^ 0.0 

Q = ] 

i7(r) = 4.9 X 10-'* 

fT(jit) — 0.060 

f7(ns) — 2.1 X 10“^ 
ujas) = 2.1 X 10“^ 

= 1.9 X 10-** 

0 

(7(A/^) = 44 
o-(Veff) - 0.029 
a(w) — 0.11 
ct(ii;o) - 0.27 
a{wa) — 0.58 

O' = 0 

a(r) = 1.2 X 10-"* 
a{nt) = 0.027 
alnf — 2.1 X 10-^ 
<t(Qs) = 2.1 X 10-** 
(7(f!|,) = 1.9 X 10-** 

16 

— 44 

o-(Veff) - 0.029 
a('w) — 0.11 
cT(4t;o) - 0.27 
a{wa) — 0.58 

a = 0 

^(r) = 2.2 X 10-'* 
(T(jit) — 0.040 
crlnf) = 2.2 x 10-** 
<T(aa) = 2.2 X 10-** 

= 2.0 X 10-** 

40 

a{AIi,) — 46 
(T{Nes) - 0.029 
a(w) - 0.12 
(j{wq) — 0.26 
(T(l[Ja) — 0.54 

Q = 0 

(7(r) - 2.1 X 10“^ 

— 0.040 

f7(71s) — 2.1 X 10“^ 
f7(Qs) - 1.9 X 10“^ 
fj(Qk) - 1.9 X 10“^ 

38 

a{AI^) — 36 
a{N^s) - 0.025 
fT(Tr;) — 0.098 

ct{wo) — 0.11 
a{wa) = 0.16 

sync-l-dust 

A - 2.5 

r-eff - 1.1 X 10“^ 

Q — 1 

a{r) = 5.2 x IQ-'* 

fT(jit) — 0.068 

alnf) = 2.2 x 10“** 
t7(as) = 2.2 X 10“** 
i7(f!|,) = 2.1 X 10-** 

0 

a{AU) - 48 
o-(Veff) - 0.032 
a(w) — 0.12 
cT(it;o) — 0.33 
a{wa) — 0.71 

a = 0 

a{r) = 1.8 X 10-“* 
a{nt) = 0.037 
alnf) = 2.2 X 10-** 
<t(Qs) = 2.2 X 10-** 
(7(f!|,) = 2.1 X 10-** 

24 

a{AU) — 48 
o-(Veff) - 0.032 
a(w) — 0.12 
c7(t^’o) — 0.33 
f7(lPa) — 0.71 

a — 0 

^(r) = 2.5 X 10-'* 
(T(jit) — 0.046 
(t(iIs) = 2.2 X 10-** 
(T(as) = 2.3 X 10-** 

= 2.1 X 10-** 

40 

(t{AU) - 49 
(T{Nes) - 0.032 
(7(ir;) — 0.13 
(j{wq) — 0.29 
f7(ltJa) — 0.59 

a — 0 

(7{r) - 2.4 X 10“^ 
(T(nt) — 0.045 
(7(71s) — 2.1 X 10“^ 
f7(Qs) - 2.0 X 10“^ 
fj(Qk) - 2.0 X 10“^ 

38 

a{AU) - 38 
f7(Veff) - 0.027 
a{w) = 0.11 
fT(itJo) — 0.12 

O'(lCa) — 0.17 

sync-t-dust 
A-expansion 
+ C-BASS 

A - 25 

r-eff — 1.8 X 10“^ 

Q — 1 

i7(r) = 9.3 X 10-'* 
fT(jit) — 0.092 
t7(lla) = 2.2 X 10-** 
t7(as) = 2.7 X 10-** 
i7(fi|,) = 2.5 X 10-** 

0 

a{AU) - 59 
o-(Veff) - 0.048 
ct(ii;) — 0.15 
cT(it;o) = 0.56 
a{wa) = 1.1 

0 = 0 

a{r) = 7.4 x 10-'* 

— 0.081 
= 2.2 X 10-** 
<t(Qs) = 2.7 X 10-** 
(7(f!k) = 2.5 X 10-** 

58 

a{AI,y) — 59 
o-(Veff) - 0.048 
cT(Tr;) — 0.15 
c7(ifo) — 0.56 
a{wa) = 1.1 

Q = 0 

(T(r) = 6.6 X 10-'* 
<j{nt) = 0.077 
alnf) = 2.2 x 10-** 
<T(aa) = 2.7 X 10-** 

= 2.5 X 10-** 

41 

a{A'U) - 57 
^{^es) — 0.045 
(t{w) = 0.16 
(j{wq) — 0.38 
C7(lt)a) = 0.71 

a — 0 

i7(r) = 6.5 X 10-'* 

= 0.077 

a{nf) = 2.2 x IQ-** 
<t(Qs) = 2.2 X 10-** 
(T(fit) = 2.2 X 10-** 

39 

a{AU) - 42 
f7(Vefr) - 0.036 
a{w) — 0.12 
ct{wo) — 0.13 
= 0.20 
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LiteBIRD-ext 

Delensing option —> 

4- comp. sep. option 

no 

delensing 

CMB X 

CMB 

CMB X 

CIB 

CMB X 

LSS 

cmb-only 

A = 1.0 

T-eff = 0.0 

a = 

a(r) = 3.0 X lO-"* 

= 0.052 

a{n,) = 3.7 x 10“^ 
(T(as) = 5.6 X 10“^ 
o-(fit) = 2.8 X 10-3 

1.0 

a{AU) = 72 
= 0.15 
a{w) = 0.16 
cr(4Co) = 0.95 
(T{Wa) = 1.9 

O' = 

tT{r) = 2.7 X 10-3 
D-{nt) = 0.049 
a(nj) = 3.7 x 10-3 
cr{Qs) = 5.6 X 10-3 
alflt) = 2.8 X 10-3 

0.79 

o-(A/^) = 72 
fr(Aeff) = 0.15 
(r(w) = 0.16 
(t(i(;o) = 0.95 

O-jlCa) = 1.9 

a = 0.-^ 

a{r) = 2.1 X 10-^ 
(T(nt) = 0.044 
(7(ns) = 3.4 X 10“^ 
cr(Q;s) = 4.6 X 10“^ 
(7(V2k) = 2.5 X 10"^ 

15 

cr{Mj.) = 61 
a(Neff) = 0.11 

a(w) = 0.15 
cr{ino) = 0.39 
ff{ina) = 0.76 

a = 0 

o-(r) = 2.1 X 10-3 
= 0.045 

tr(ns) = 2.9 X 10-3 
tr(as) = 2.3 x 10-3 
tr(rit) = 2.3 X 10-3 

46 

a(AU) = 45 
(T(Aeff) = 0.074 
(t(w) = 0.11 
ai'ico) = 0.12 
a’(wa) = 0.21 

sync+dust 

A = 1.7 

reff = 1.0 X 10-^ 

a = 

a(r) = 5.8 x lO-"* 

= 0.066 
a(n,) = 3.9 x 
fT(as) = 5.8 X 10-3 
o-(S4) = 3.1 X 10-3 

1.0 

a{AU) = 80 

= 0.16 
(T(itf) = 0.17 
crfico) = 1.0 

= 2.1 

O' = 

a(r) = 5.2 x 10-3 
D-{nt) = 0.064 
cr(ns) = 3.9 X 10-3 
cr{Qs) = 5.8 X 10-3 
a-{nt) = 3.1 X 10-3 

0.83 

fr(A4) = 80 
fr(Aeff) = 0.16 

cr(w) = 0.17 
cr(iuo) = I'O 

0-(l(;g) = 2.1 

a = Os 

a{r) = 3.8 x lO”'* 
(T(nt) = 0.057 
cr{ns) = 3.5 x 10“^ 
cr(Q;s) = 4.7 X 10“^ 
(7(Qk) = 2.6 X 10-^ 

17 

a(M^) = 62 
ff(Neff') = 0.11 

(7{it;) = 0.15 
ff{ino) = 0.39 
ar(Wa) = 0.76 

a = 0 

tr(r) = 3.8 x 10-3 
'T(nt) = 0.057 
(r(ns) = 3.0 X 10-3 
tT(as) = 2.3 X 10-3 
cr(nt) = 2.3 X 10-3 

47 

(t(A 4) = 46 
(T(Aeff) = 0.076 
(t(w) = 0.11 
(t(ico) = 0.13 
a’(wa) = 0.21 

sync+dust 
A-expansion 
+ C-BASS 

A = 7.8 

T-eff = 6.3 X 10-^ 

Q = 

a(r) = 1.2 X 10-3 
= 0.087 

o-(ns) = 4.9 X 10-3 
a(a„) = 6.8 x 10-3 
(T(nt) = 5.5 X 10-3 

1.0 

cr(M„) = 1.4 X 103 
cr(Ar,,a.) = 0.20 

(t{w) = 0.29 
frjico) = 1.3 
crlwa) = 2.7 

O' = 

(7(r) = 1.2 X 10“^ 
(7{nt) = 0.087 
cr{ns) = 4.9 x 10"^ 

cr{as) = 6.8 X 10“^ 

(7{nk) = 5.5 X 10“^ 

0.94 

tT(M^) = 1.4 X 10^ 

fr(Aeff) = 0.20 

cr(w) = 0.29 
o-(wo) = 1.3 
(T(Wa) = 2.7 

Q = 0. 

a{r) = 1.0 X 10“^ 
cr(nt) = 0.084 
(7(ns) = 4.1 X 10"^ 
cr(Qs) = 5.0 X 10“^ 
crfOk) = 2.7 X 10“^ 

>1 

cr(AU) = 66 

0-{Aeff) = 0.12 

cr(w) = 0.16 
cr(wo) = 0.40 
cr(w^) = 0.78 

a = 0 

tr(r) = 1.0 X 10-3 
= 0.084 

cr(iis) = 3.3 X 10-3 
ar(as) = 2.3 x 10-3 
tr(fik) = 2.5 X 10-3 

52 

(T(Af^) = 49 

a(N^ff) = 0.086 
a(w) = 0.12 
cr(ico) = 0.13 

cr(n;a) = 0.22 

sync+dust 
A-expansion 
+ QUIJOTE 

A = 7.3 

T-eff = 5.1 X 10-^ 

a = 

a(r) = 1.1 X 10-3 
= 0.086 

o-(na) = 4.9 X 10-3 
fT(as) = 6.7 X 10-3 
= 5.3 X 10-3 

1.0 

a(M„) = 1.4 X 103 
cr(Afeft) = 0.20 
(t{w) = 0.28 
crjico) = 1.3 

<7(tCa) = 2.6 

O' = 

CT{r) = 1.1 X 10-3 
cr(nt) = 0.086 
a-{ns) = 4.9 x 10-3 
cr(as) = 6.7 X 10-3 
= 5.3 X 10-3 

0.94 

tx(M^) = 1.4 X 10^ 

O-(Aeff) = 0.20 

a(w) = 0.28 

o-jico) = 1'3 

cr(Wa) = 2.6 

Q = 0. 

a{r) = 9.6 x 10-3 
= 0.082 

cr(ns) = 4.1 X 10-3 
cr(Qs) = 5.0 X 10-3 
crfslt) = 2.7 X 10-3 

>1 

cr(AU) = 66 
^{A^eff) = 0.12 
(t(w) = 0.16 
cr(wo) = 0.40 
CT(Wa} = 0.78 

a = 0 

tr(r) = 9.6 X 10-3 
= 0.082 

a-(ns) = 3.3 x 10-3 
o-(as) = 2.3 X 10-3 
ir(fik) = 2.5 X 10-3 

52 

cr(Af^) = 48 

(T(Aeff) = 0.086 

a(w) = 0.12 
(t(ico) = 0.13 

(T(Wa) = 0.22 


COrE+ 

Delensing option —> 
i comp. sep. option 

no ]| CMB X 

delensing jj CMB 

CMB X 

CIB 

CMB X 

LSS 

cmb-only 

A - 1.0 

T-eff ^ 0.0 

a = ] 

u(r) = 2.9 X 10-3 
f7(7it) — 0.048 

I7(lla) = 1.9 X 10-3 
u(as) = 2.1 X 10-3 
u(fik) = 1.7 X 10-3 

0 

f7(Af^) - 38 
(T(N^ff) - 0.036 

(7(w) = 0.100 

fT(i[Jo) — 0.39 
^(iPa) — 0.79 

a = 0 

u(r) = 1.7 X 10-3 
(T(jit) — 0.034 

= 1.9 X 10-3 
(jJqs) = 2.1 X 10-3 
u(Qt) = 1.7 X 10-3 

33 

f7(A4) - 38 

- 0.036 
a('w) = 0.100 
(T(nJo) — 0.39 
(7(Wa) — 0.79 

a = 0 

(j(r) = 1.8 X 10-3 
cj(nt] — 0.036 
<T(ras) = 1.9 X 10-3 
(T(as) = 2.2 X 10-3 
(T(fik) = 1.7 X 10-3 

40 

(t(A4) - 39 
(7(N^ff) - 0.035 

(7(w) - 0.12 

(t(i[Jo) — 0.31 
(T(l[Ja) — 0.57 

a = 0 

a(r) = 1.8 X 10-3 
a{n,) — 0.036 
<T(ns) = 1.8 X 10-3 
a{a,) = 1.8 X 10-3 
(7(fik) = 1.6 X 10-3 

38 

(7(A4) - 32 
t7(N^ff) - 0.029 
c7(w) = 0.094 

(T(n;o) — 0.10 
o-(wa) — 0.16 

sync+dust 

A - 2.3 

T-eff - 1.4 X 10“^ 

a — ] 

<T(r) = 5.5 X 10-3 
— 0.060 

i7(raa) = 1.9 X 10-3 
i7(as) = 2.3 X 10-3 
t7(fik) = 1.7 X 10-3 

0 

fT(A/^) - 40 
a(Neff) - 0.042 
cr(w) = 0.10 
fT(i[Jo) — 0.45 
^(iPa) — 0.92 

a = 0 

cr(r) = 3.4 X 10-3 
(T(jit) — 0.048 

= 1.9 X 10-3 
(jJqs) = 2.3 X 10-3 
(7(f!k) = 1.7 X 10-3 

45 

fT(A4) - 40 
CT(AefF) - 0.042 
a(w) = 0.10 
(t(i(Jo) — 0.45 
= 0.92 

a = 0 

(T(r) = 3.3 X 10-3 
= 0.047 
= 1.9 X 10-3 
(T(as) = 2.3 X 10-3 
(T(S2k) = 1.7 X 10-3 

40 

(T(Afi,) — 41 
a(Neff) - 0.039 
a(w) - 0.12 
(t(i[io) — 0.32 
— 0.60 

a = 0 

a{r) = 3.1 X 10-3 
fT(7it) — 0.046 
r7{nj) = 1.9 X 10-3 
a{as) = 1.8 X 10-3 
(jffik) = 1.6 X 10-3 

38 

(7(A4) - 33 
a(Neff) - 0.031 
cr(w) — 0.096 

(T(n;o) — 0.11 

— 0.17 

sync+dust 
A-expansion 
+ C-BASS 

A - 9.8 

T-eff — 3.9 X 10“® 

a — ] 

u(r) = 1.0 X 10-3 
f7(7it) — 0.076 

I7(lla) = 1.9 X 10-3 
u(as) = 2.6 X 10-3 
u(fik) = 2.0 X 10-3 

0 

a(AU) - 47 
a(Neff) - 0.054 
a(w) — 0.12 
cT(it)o) = 0.58 
alwa) — 1.2 

a = 0 

u(r} = 8.9 X 10-3 
{T(nt) — 0.073 
a(ns) = 1.9 X 10-3 
cr(aj) = 2.6 X 10-3 
(7(f!k) = 2.0 X 10-3 

68 

f7(A4) - 47 
o-(Aeff) - 0.054 
( t ( w ) — 0.12 
(t(i[;o) — 0.58 

(7(u!a) — 1.2 

a = 0 

(T(r) = 7.8 X 10-3 
(T(iit) = 0.070 
<T(ras) = 1.9 X 10-3 
(T(as) = 2.6 X 10-3 
(T(S2k) = 1.8 X 10-3 

41 

(7(Af,Q - 43 
o-(Aeff) - 0.049 
a(w) — 0.13 
(t(i[Jo) — 0.34 
(T(l[Ja) — 0.64 

a = 0 

<j(r) = 7.7 X 10-3 
a{nt) — 0.069 
<T(ns) = 1.9 X 10-3 
a(a,) = 1.9 X 10-3 
(7(fik) = 1.7 X 10-3 

39 

a-(A/i,) — 35 
a(N^ff) - 0.037 
a-(w) — 0.099 
(T(njo) — 0.11 
— 0.18 

sync+dust 
A-expansion 
+ QUIJOTE 

A - 7.4 

T-eff - 1.9 X 10“*= 

a — ] 

t7(r) = 8.5 X 10-3 
o-(raj) = 0.072 
fT(jls) — 1.9 X 10-3 
i7(as) = 2.5 X 10-3 
t7(f!k) = 1.9 X 10-3 

0 

ct(A/^) — 45 
cr(Neff) - 0.051 
a(w) — 0.12 
cT(it;o) = 0.55 

(7(Wa) = 1.1 

a = 0 

<T(r) = 7.2 X 10-3 
cr(nt) = 0.067 

— 1'9 X 10-3 
a(as) = 2.5 x 10-3 
<T(rik) = 1-9 X 10-3 

64 

a(AI^) — 45 
a(Neff) - 0.051 
o-(w) — 0.12 
cT(if)o) = 0.55 
a(wa) = 1.1 

a = 0 

(7(r) - 6.3 X 10“^ 
(T(nt) — 0.064 
(7(jis) — 1.9 X 10“^ 
(7(as) = 2.5 X 10“^ 
(jff^k) - 1.8 X 10“^ 

41 

(7(A/,Q - 42 
(7(Aeff) - 0.047 
(7(w) — 0.13 
(t(i[io) — 0.34 
cT(it;a) = 0.63 

a = 0 

a{r) = 6.2 X 10-3 
fT(7it) — 0.064 
fj(ns) — 1.9 X 10-3 
(7(as) = 1.9 X 10-3 
<T(rit) = 1.7 X 10-3 

39 

a(A/i,) — 34 

- 0.036 
— 0.098 

(T(n;o) — 0.11 
C7(Wa) = 0.18 
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PIXIE 


Delensing option — ^ 

4- comp. sep. option 

no 

delensing 

CMB X 

CMB 

CMB X 

CIB 

CMB X 

LSS 

cmb-only 

A - 1.0 

Peff ^ 0.0 

a = 

a(r) = 3.1 X 10-* 
o-(n,) = 0.070 
a(nj} = 0.070 
<t{Qs) = 0.034 
o-(r!t) = 0.018 

- 1.0 

a{AU) - 4.6 X 10^ 
a(7Veff) - 1.3 
cr{w) = 0.61 
(^(ico) = 2.5 
(T{Wa) — 5.6 

a = 

o-(r) = 3.0 X 10“.* 
(j{nt) — 0.070 
cr(nf) = 0.070 
cr(as) = 0.034 
(T(rik) = 0.018 

1.00 

a(AU) - 4.6 X 10^ 

^ 1-3 

a(w) = 0.61 
c7(wo) = 2.5 
— 5.6 

a = 0. 

cT(r) - 2.9 X 10^^ 
fT(nt) — 0.069 
a-(ns) — 6.3 x 10“^ 
ct(qs) - 7.6 X 10“^ 
o-(nk) - 3.5 X 10“^ 

15 

(7(A'U) ^ 86 

^ 0.72 
cr(w) = 0.18 
fr(i[Jo) — 0.63 
^ 1,3 

£1 = 0. 

£7(r) = 2.9 X 10“'* 
£T(nt) = 0.067 
£7+) = 4.8 X 10“** 
£7+) = 4.0 X 10“** 
£7(Slk) = 3.3 X 10“** 

34 

a(A4) - 64 

a(N^ff) - 0.68 

(t(w) = 0.14 
(t(ico) — 0.21 
(T(Wa) = 0.40 

sync+dust 

A - 1.2 

reff = 1.4 X 10“^ 

a = 

a(r) = 5.5 x 10“"* 
(T(nt) — 0.090 
= 0.070 
o-{q,) = 0.034 
<T(r!k) = 0.018 

- 1.0 

(t{AU) - 4.7 X 10^ 

^ 1.3 
a-{w) = 0.63 
^■(ico) = 2.5 
(T(tCa) = 5.7 

a = 

a{r) = 5.4 x 10“.* 
{T{nt) — 0.090 
cr(nf) = 0.070 
cr(as) = 0.034 
tT(rit) = 0.018 

1.00 

(7(AU) - 4.7 X 10^ 
t7(iVeff) - 1.3 
(7(w) = 0.63 
cr(wo) = 2.5 
— 5.7 

a-0. 

o-(r) - 5.1 X 10“^ 
CT(rit) — 0.089 
f7(ns) — 6.4 X 10“^ 
o-(qs) - 7.6 X 10“^ 

- 3.6 X 10“3 

35 

o-(A4) - 87 
cT(N^ff} - 0.74 

(7(w) — 0.18 

(7(iuo) — 0.63 
f7(l[Ja) — 1.3 

a - 0. 

(7(r) - 5.1 X 10-^ 
f7(nt} — 0.088 
c7(ns) — 4.9 X 10-3 
cT(as) - 4.1 X 10-3 
f7(flk) - 3.3 X 10-3 

34 

a-(AIi,) — 65 
(T(Aeff) — 0.69 
(t(w) = 0.14 
o-(wo) — 0.21 
(7(Wa) = 0.41 

sync+dust 
A-expansion 
+ C-BASS 

A - 1.6 

J-eff — 1.0 X 10 ^ 

a = 

(T{r) = 5.9 X 10-* 
a-(nt) — 0.092 
<T(tis) = 0.070 
<t(q,) = 0.034 
(7+) = 0.019 

- 1.0 

a{AU) - 4.8 X 10^ 
o-iWeff) ^ 1.3 
a(w) = 0.65 
(T(4t;o) = 2.5 

CT{u'a) — 5.8 

Q = 

tT(r) = 5.9 X 10“.* 
cj(nt) — 0.092 
cr(nf) = 0.070 
(7(as) = 0.034 
(T(rik) = 0.019 

1.00 

a(AI^) - 4.8 X 10^ 

^ 1.3 
ct(w) = 0.65 
cT(ido) = 2.5 
— 5.8 

a-0. 

a(r) - 5.6 x 10^^ 
fT(nt) — 0.091 

CT(ns) — 6.4 X 10“^ 
(T(as) - 7.7 X 10'^ 
o-(nk) - 3.6 X 10-3 

35 

a(AU) - 88 
o-(Aeff) - 0.75 
— 0.18 
(7(i[Jo) — 0.64 
(7(Wa) — 1.3 

Q - 0. 

a(r) - 5.5 x 10-^ 
(T(nt) — 0.090 
(T(ns) — 4.9 X 10-3 
(t(q's) — 4.1 X 10-3 
fT(flk) - 3.4 X 10-3 

34 

(T(A'fi,) — 65 
a(N,ff) - 0.70 
cr(w) = 0.14 
(7(wo) = 0.22 
(T(Wa) = 0.41 

sync+dust 
A-expansion 
+ QUIJOTE 

A - 1.4 

T-eff - 9.7 X 10'^ 

a = 

o-(r) = 5.7 X 10-* 
(T(nt) — 0.091 
o-(n,) = 0.070 
o-{q,) = 0.034 
<T(r!t) = 0.019 

- 1.0 

^ 4.8 X 10^ 
'T(Aeff) ^ 1.3 
(t{w) = 0.64 
£7(4Co) = 2.5 
(T(4Ca) = 5.7 

a = 

(T(r) = 5.7 X 10“.* 
{T(nt) — 0.091 
cr(nf) = 0.070 
(T(as) = 0.034 
tr(rik) = 0.019 

1.00 

(7(M^) - 4.8 X 10^ 
t7(iVeff) - 1.3 
cr(w) = 0.64 
cr(wo) = 2.5 
= 5.7 

a-0. 

a(r) - 5.4 x 10“^ 
a(rit) — 0.090 
a(ns) = 6.4 x IQ-^ 
a(as) - 7.7 x lO'^ 
c7(flk) = 3.6 X 10-3 

35 

(7(AU) - 88 
o-(Aeff) - 0.75 
(7(w) — 0.18 

(7(iuo) — 0.63 
a(wa) = 1.3 

a - 0. 

c7(r) - 5.3 X 10-^ 
f7(nt) — 0.089 
cT(ns) — 4.9 X 10-3 
cT(as) - 4.1 X 10-3 
cT(i7k) = 3.3 X 10-3 

34 

— 65 

cr(Neff) - 0.70 
(t(u;) = 0.14 
o-(wo) — 0.22 
(7(Wa) = 0.41 


LiteBIRD-ext x COrE+ 

Delensing option —> 

1 comp. sep. option 

no 

delensing 

CMB X 

CMB 

CMB X 

CIB 

CMB X 

LSS 

cmb-only 

A - 1.0 

Peff ^ 0.0 

Q = ] 

£7(r) = 2.8 X 10“'* 
£7(iit) = 0.047 
£7(jla) = 1.9 X 10“** 
£7(as) = 2.1 X 10“** 
£7(f!k) = 1.7 X 10“** 

0 

fT(A4) - 38 
o-(A^eff) - 0.036 

(7(w) = 0.100 

cT(it;o) — 0.39 
fT(TUa) — 0.79 

a = 0 

£7(r) = 1.5 X 10“'* 
£7(71t) = 0.032 
£7(jls) = 1.9 X 10“** 
£7 (Qs) = 2.1 X 10“** 
£7(r!k) = 1-7 X 10“** 

33 

a(AU) - 38 
a(N^ff) - 0.036 
a(w) = 0.100 
o'(T.fo) — 0.39 
o'(iPci) — 0.79 

a = 0 

£7(r) = 1.7 X 10“'* 
£7(ra,) = 0.035 
£7(ras) = 1.9 X 10“** 
£7(aa) = 2.2 X 10“** 
£7(f!k) = 1-7 X 10“** 

40 

(t(A 4) - 39 
f^(A'eff) — 0.035 
(7(w) - 0.12 
(7(wo) — 0.31 
(T(l[Ja) — 0.57 

a = 0 

£7(r) = 1.6 X 10“'* 
£7(nt) = 0.035 
£7(71s) = 1.8 X 10“** 
£7(Qs) = 1.8 X 10“** 
£7(fit) = 1-6 X 10“** 

38 

a{AU) - 32 
(7(Aeff) - 0.029 
a{w) = 0.094 

ct{wo) — 0.10 
a{wa) = 0.16 

sync+dust 

A - 1.8 

r-eff - 3.5 X 10-3 

Q — ] 

£7(r) = 4.7 X 10“'* 
£7(jit) = 0.057 
£7(raa) = 1.9 X 10“** 
£7(as) = 2.2 X 10“** 
£7(f!k) = 1.7 X 10“** 

0 

fT(A/^) - 40 
fT(AefF) - 0.040 
(7(w) = 0.10 
cT(it;o) — 0.43 

(T(Wa) — 0.88 

a = 0 

£7(r) = 2.5 X 10“'* 
£7(71t) = 0.042 
£7(jls) = 1.9 X 10“** 
£7(Qs) = 2.2 X 10“** 
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